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Abstract 



We establish the self-averaging properties of the Wigner transform of a mixture of states in 
the regime when the correlation length of the random medium is much longer than the wave 
length but much shorter than the propagation distance. The main ingredients in the proof are 
I the error estimates for the semiclassical approximation of the Wigner transform by the solution 

. of the Liouville equations, and the limit theorem for two-particle motion along the characteristics 

of the Liouville equations. The results are applied to a mathematical model of the time-reversal 
experiments for the acoustic waves, and self-averaging properties of the re-transmitted wave are 
proved. 



1 Introduction 

1.1 The Wigner transform of mixtures of states 



O 

' The Wigner transform is a useful tool in the analysis of the semi-classical limits of non-dissipative 

evolution equations as well as in the high frequency wave propagation |jl^, |l^, |2^. It is defined 
! as follows: given a family of functions /e(t,x) uniformly bounded in L°°([0, T]; L^(M'^)) its Wigner 

transform is 

^ TV.(t,x,k)= |e^k.y^,(t,x-^)/;(i,x+^)^. (1.1) 

^ . The family is uniformly bounded in the space of Schwartz distributions 5'(M'^ x M*^), and all its 

X 



limit points are non-negative measures of bounded total mass [14, 17|. It is customary to interpret a 
■ limit Wigner measure W as the energy density in the phase space, since the limit points of = 

are of the form n{t,x) = J W{t,x,'k)dk provided that the family is e-oscillatory and compact at 
infinity [|l4|. However, while neither nor its limit n(t,x) satisfy a closed equation, both Ws and W 
usually obey an evolution equation when the family /(.(tjx) arises from a time-dependent PDE. This 
makes the Wigner transform a useful tool in the study of semiclassical and high frequency limits, 
especially in random media |[l|, |2|, |10|, 18, |2^. However, a priori bounds on the Wigner transform 
other than those mentioned above are usually difficult to obtain. It has been observed in that 
the Wigner transform of a mixture of states 



£ 



Wsix,k)= |e^k.y^,(x-^;C)/;(x+^;C) 
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enjoys better regularity properties. The family above depends on an additional "state" parameter 
C ^ S, where 5 is a state space equipped with a non-negative bounded measure dfi{(). Typically this 
corresponds to introducing random initial data for at t = and estimating the expectation of 



with respect to this randomness. This improved regularity has been used, for instance, in pq , 22] 
in the analysis of the average of the Wigner transform of mixtures of states in random media and 
in in order to obtain an asymptotic expansion for the Wigner transform of a mixture of states. 

The purpose of this paper is to analyze the self-averaging properties of moments of the mixed 
Wigner transform of the form J Ws{t,x,\<.)S(k)dk, where 5'(k) is a test function, and the family 
/e(t,x;(^) satisfies the acoustic equations. This problem arises naturally in the mathematical study 
of the experiments in time-reversal of acoustic waves that we will describe in detail below. How- 
ever, apart from the time-reversal application, the statistical stability of such moments provides 
an important key to understanding the physical applicability of the limit equations for the Wigner 
transform in random media in the situations when results for each realization are more relevant than 
the statistically averaged quantities. 

We start with the wave equation in dimension d > 3 

' A^ = (1.2) 



C2(x) 9t2 



and assume that the wave speed has the form c(x) = cq + \/5ci(x). Here cq > is the constant 
sound speed of the uniform background medium, while the small parameter 5 <C 1 measures the 
strength of the mean zero random perturbation ci. Rescaling the spatial and temporal variables 



X = x'/(5 and t = t' /6 we obtain (after dropping the primes) equation (1.2) with rapidly fluctuating 
wave speed 

C5(x) = CO + V6ci (I) . (1.3) 

It is convenient to re- write (11.21) as the system of acoustic equations for the "pressure" p = —(f>t and 

— c 
"acoustic velocity" u = —Vcj): 

^ + V(c,(x)p) = (1.4) 
^ + C5(x)V-u = 0. 

The energy density for ( |1.4D is £'(t,x) = |up +p'^: J E{t,x)dx = const is independent of time. We 
will denote for brevity v = (u,p) G C^"^^ and write (|l.4| ) in the more general form of a first order 
linear symmetric hyperbolic system 

^ + (^5 (x)vf (x)) = 0. (1.5) 

In the present case, the symmetric matrices As and are defined by 

^^(x) = diag(l, 1, l,C5(x)), and D-' = (g) e^+i + e^+i ® e^, j = 1, . . . , d. (1.6) 

Notice that the matrices are independent of x. Here G M.'^^^ is the standard orthonormal 
basis: {em)k = ^mk- The dispersion matrix for ( |1.5| ) is 



Po^(x, k) = iAs{x)kjD^ As{x) = ics{x)kjD^ = ics{x) k e^+i + e^+i k , k = ^ kjej. (1.7) 
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The self-adjoint matrix {—iPq) has an eigenvalue Aq = of multiplicity d — 1, and two simple 
eigenvalues Aj_(x, k) = ibc5(x)|k|. The corresponding eigenvectors are 

b^= (k;^,0), m = l,...,d-l; b± = |^ A ± e,+i j , (1.8) 

where k^ € M.'^ is the orthonormal basis of vectors orthogonal to k. 

We assume that the initial data vo(x;C) = vf(0,x;(") = {—e'V4>Q,l/cs4>o) for ( |1.5D is an e- 
oscillatory and compact at infinity family of functions uniformly bounded in L^(]R°') [^J] for each 
"realization" ( of the initial data. The scale e of oscillations is much smaller than the correlation 
length 5 of the medium: e <^ 6 <^ 1. The (d + 1) x (d + 1) Wigner matrix of a mixture of solutions 
of (|l.5|) is defined by 

W';(t.x,k)= / .*.vJ(t,x-f;C)vf(,x+f;C)«. 

R<*x5 

The non-negative measure d/i has bounded total mass: Jgdfi{(^) < oo. It is well-known [|l^, 17| that 
for each fixed 5 > (and even without introduction of a mixture of states) one may pass to the limit 
e ^ and show that converges weakly in 5'(M'^ x W^) to 

W^'^(t,x,k) = M^(t,x,k)b+(k) ®b+(k) +nl(t,x,k)b_(k) 0b_(k). 

The scalar amplitudes satisfy the Liouville equations: 

^ + VkAi • Vx4 - VxAi • Vk4 = 0. (1.9) 



dt 



Furthermore, one may formally pass to the limit (5 — > in ( |1.9D and show that (see [^) Ejti^} 
converge to the solution of 



du± - d 

± CqK ■ VxUi 



dt dkm 
Here k = k/|k|, and the diffusion matrix D is given by 



(1.10) 



1 d'^R(cos'k) , 

Drnn = -J -^^^ds, (1.11) 

where -R(x) is the correlation function of ci: E{ci(y)ci(x + y)} = -R(x). 

The purpose of this paper is to make the passage to the limit e, 5 — > rigorous for a mixture 
of states (and eliminate the consecutive limits e — > then (5 — > 0) and establish the self-averaging 
properties of moments of the form s^(t,x) = J (t, x, k)S'(k)(ik, where S G L^(M'^) is a given test 
function. 

The assumption that e <^ 5 is formalized as follows. We let /C^ = {(e, 5) : S >\ lne|~^/'^+''}, 
with < fj, < 2/3 and assume that {e,5) G /C^ for some /i G (0,2/3). From now on, /x is a given 
fixed number in (0, 2/3). 
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1.2 The random medium 

We make the following assumptions on the random field ci(x). Let (r2,C,P) be a certain probability 
space, and let E denote the expectation with respect to P and || • \\p denote the respective norm 
for any p G [1, +00]. We suppose further that ^ M is a measurable, strictly stationary, 

mean-zero random field, that is pathwise C^-smooth and satisfies 

A := ess-sup |VxCi(x;u;)| < +00, i = 0, 1,...,4. (1-12) 

We assume in addition that ci is exponentially (/)-mixing. More precisely, for any i? > we let 
C]:j := 0"{ci(x) : |x| < R} and Cfj := cr{ci(x) : |x| > R}. We also define 

0(p) := sup[|P(i?) - FiB\A)\ : R>0,AGCn,BG Cfj+^J, 

for all p > 0. We suppose that there exists a constant Ci > such that 

(j){p) < le-'^'P , Vp>0. (1.13) 

We let also 

i?(y) =E[ci(y)ci(0)], y G M'^ 

be the covariance function of the field ci(-) and note that ( 1.13| ) implies that there exists a constant 
C2 > such that 

|V^i?(y)| < C72e-|yl/^2^ Vy G M*^, m = 0, • • • , 4. (1.14) 

Finally we assume that R G C°°(M'^), this condition will be used only to establish the hypoellipticity 
of ( [1.10 ). Notice that sufficiently regular random fields with finite correlation length satisfy the 



hypotheses of this section. The exponential ^-mixing assumption was used in [15| to analyze the 



solutions of Liouville equations with random coefficients. Their techniques lay at the core of our 



proof of the mixing properties presented in our main result. Theorem |1.1| , below. 
1.3 The main result 

We assume that the initial Wigner transform W^{0, x, k) is uniformly bounded in L^(]R'^ x M*^) and 
W^{0, X, k) ^ VFo(x, k) strongly in L^{R'^ x M'^) as /C^ 3 (e, 6) 0. (1.15) 
We also assume that the limit Wo G Cc(M'^ x W^) with a support that satisfies 

supp VFo(x,k) C X = {(x,k) : |x| < C, C-^<\k\<C}. (1.16) 



Note that ( p.. 151 ) may not hold for a pure state since ||W'£||2 = (27re) ''^^ll/elli 0- We will later 



present examples where it does hold for a mixture of states. Furthermore, we assume that Wq has 
the form 

VFo(x,k) =u°b+®b+ + u°b_®b_ (1.17) 

and let 

W(t,yi,k) = ii+(t,x,k)b+(k) ®b+(k) +«_(*, x,k)b_(k) ®b_(k). (1.18) 

The functions u± satisfy the Fokker-Planck equation ( |1.10D with initial data as in ( [1.171 ). The 
main result of this paper is the following theorem. 
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Theorem 1.1 Let us assume that the random field ci(x) satisfies the assumptions given in Section 
|r| and that the initial data W/(0,x,k) satisfies and Let 5(k) G L'^{M'^) be a test 

function, and define the moments 



W^{t,x,k)S{k)dk and s{t,x) 



W{t,x,k)S{k)dk, 



where W is given by { l.lt ). Then for each t > we have 

eI^J |sf(t,x) - s(t,x)p(ix| ^ 



(1.19) 



as K,^ 3 



Theorem |1.1| means that the moments converges to a deterministic hmit. The main apphcation 



of Theorem 1.1 we have in mind is the mathematical modehng of refocusing in the time-reversal 
experiments we present in Section ^. 

Our results may be generalized in a fairly straightforward manner to other wave equations that 
may be written in the form (|l.5| ), which include acoustic equations with variable density and com- 



pressibility, electromagnetic and elastic equations [20|. 

The papers is organized as follows. The mathematical framework of the time-reversal experi- 
ment as well as the main result concerning the self-averaging properties of the time reversed signal, 
Theorem 2.1, are presented in Section |2[ Section ^ contains the derivation of the Liouville equations 
in the L^-framework. Some straightforward but tedious calculations from this Section are presented 
in Appendices ^ and The limit theorem for the two-point motion along the characteristics 
of the Liouville equations, Theorem 4^, is presented in Section ^ Theorem LI follows from this 
result. The proof of Theorem 4.4 is contained in Section ^. 
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2 Refocusing in the time-reversal experiments 

2.1 Mathematical formulation of the time-reversal experiments 

Refocusing of time-reversed acoustic waves is a remarkable mathematical property of wave prop- 
agation in complex media that has been discovered and intensively studied by experimentalists in 
the last decade (see ||T^, |l6| and also [^] for further references to the physical literature). A typical 
experiment may be described schematically as follows. A point source emits a localized signal. The 
signal is recorded in time by an array of receivers. It is then reemitted into the medium reversed 
in time so that the part of the signal recorded last is reemitted first and vice versa. There are two 
striking experimental observations. First the repropagated signal tightly refocuses at the location 
of the original source when the medium is sufficiently heterogeneous even with a recording array 
of small size. This is to be compared with the extremely poor refocusing that would occur if the 
heterogeneous medium were replaced by a homogeneous medium. Second, the repropagated signal 
is self-averaging. This means that the refocused signal is essentially independent of the realization 
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of a random medium with given statistics, assuming that we model the heterogeneous medium as a 
random medium. 

The first mathematical study of a time-reversal experiment has been performed in in the 



framework of one-dimensional layered random media. The one-dimensional case has been further 
studied in [^], and a three-dimensional layered medium was considered in |13]. The time-reversal 



experiments in an ergodic domain have been analyzed mathematically in The basic ideas that 
explain the role of randomness in the refocusing beyond the one-dimensional case were first outlined 
in in the parabolic approximation of the wave equation, that was further analyzed in H, 19]. 



Time-reversal in the general framework of multidimensional wave equations in random media has 
been studied formally in |^ Q. One of the purposes of this paper is to present the rigorous proof of 
some of the results announced in Q. 

The re-transmission scheme introduced in |^, ^ is as follows. Consider the system of acoustic 
equations (|1.4D (or, equivalently, (L5) for the pressure p and the acoustic velocity u(t,x)). The 



initial data for ( |1.5D is assumed to be localized in space: 



Here xq G M is the location of the original source, and Sq G 5(M ) is the source shape function. 
The small parameter e <C 1 measures the spatial localization of the source. The signal v^(t, T) is 
recorded at some time t = T, processed at the recording array and re-emitted into the medium. The 
new signal is the solution of (|1.5D on the time interval T <t < 2T with the Cauchy data 



vf(r,x) = r[/,*[xvf(r)](xMx). (2.2) 

The initial data ( |2.2| ) reflects the process of recording of the signal at the array and its smoothing by 
the recording process. The kernel /e(x) = e~^/(x/e) represents the smoothing. The array function 
x(x) is either the characteristic function of the set of the receivers, or some non-uniform function 
supported on this set. We will assume for simplicity that /(|y|) is radially symmetric, and, moreover, 

X e Cc(M"'), / G CciM.'^), supp/(k) C {0 < C-^ < |k| < C < 00} (2.3) 

where 

/(k) = I e~'^-yf{y)dy 

is the Fourier transform of /. The matrix T corresponds to the linear transformation of the signal. 
The pure time-reversal corresponds to keeping pressure unchanged but reversing the acoustic velocity 
so that r = To := diag(— 1, —1, —1, 1). However, this is only one possible transformation, and while 
we restrict T to the above choice our results may be extended to more general matrices F, or even 
allow F be a pseudo-differential operator of the form F(x, eD). 

The re-propagated field near the source at time t = 2T is defined as a function of the local 
coordinate ^ and of the source location xq : 

vP(^;xo)=vf(2T,xo+eO- 

2.2 The re-propagated signal and the Wigner transform 

Let us assume that the random field ci satisfies the assumptions of Theorem outlined in Section 



1.2. Then Theorem 1.1 implies the following result. 
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Theorem 2.1 Under the assumptions made above the re-propagated field vf'^(^,xo) converges as 
/C^ 9 (e, <5) ^ 

dk 



v^(^,xo) = I e^i'-^[u+(r,xo,k)(5o(k),b_(k))b+(k)+n_(T,xo,k)(5o(k),b+(k))b„(k)] 



(2vr) 



d 



in the sense that 



supe|/ |vP(^,xo)-v^(txo)pdxo| ^0. (2.4) 

The functions u±{t,x,'k) are the solutions of the Fokker-Planck equation Ijil.K] ) with initial data 
a±(0,x,k) = |x(x)|2/(k). 



The proof of Theorem |2.1| is based on Theorem 1.1 and a representation of the re-propagated signal 



in terms of the Wigner transform of a mixture of solutions of the acoustic wave equations. The 
latter arises as follows. Let Qf(i,x;q) be the matrix- valued solution of (pT^) with initial data 



Qf(0,x;q)=x(x)e*'i-/^/, (2.5) 

where / is the (d + 1) x (d + 1) identity matrix, x(x) is the array function, and q S is a fixed 
vector. It plays the role of the "state" of the initial data. Physically Qf describes evolution of a 
wave that is emitted by the recorders-transducers with a wave vector q. The Wigner transform of 
the family Qf(i,x;q) is 

M^,^(t,x,k;q) = J e'k-yQ5(t,x-|^;q)Qf (t,x+^;q)^. (2.6) 

The corresponding "mixed" Wigner transform is 

W^{t, X, k) = I W,'{t, X, k; q)/(q)dq. (2.7) 

Then the re-propagated signal is described as follows in terms of W^. 
Lemma 2.2 The re-propagated signal may be expressed as 

vf'^(l,xo) = I e^^<^-y^W,\T,xo + ^^,k)roSo(y)^. (2.8) 



Proof. Let G{t, x; y) be the Green's matrix of ( |1.5D , that is, solution of ( |1.5D with the initial data 
G(0, x;y) = 16{x — y). Then the signal arriving to the recorders-transducers array is 

vf(r,x) = J G(r,x;y)vf(0,y)dy = J G(r,x;y)So (^^^) dy 



and the re-emitted signal is 



vf(T,z)= J ro/e(z-z')x(z)x(z)vf(r,z')dz'. 



Therefore we obtain 



vf'^(^,xo) = J G{T,xo + e$-z)^,{T,z)dz (2.9) 
= / G(T,xo +e^;z)ro/ (^) x(z)x(z')G(T, z'; y)So [^'^ '^'^''^ 
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However, we also have 

roG(t,x;y)ro = G*(t,y;x). 
This is seen as foUows: a solution of (|1.5[ ) satisfies 



(2.10) 



v(i,x) = y" G{t - s,x;y)v{s,y)dy 



for all < s < t. Differentiating the above equation with respect to s, using ( |1.5D for v(s, y) and 
integrating by parts we obtain 







■ ^^^^ g"'""'^^ + ^ m - s, x; y)Asiy)) D^As{y) ] v(s, y). 



Passing to the limit s ^ and using the fact that the initial data vo(y) is arbitrary we obtain 

(2.11) 



- ^ (G(t,y;x)^5(x))Z)^A,(x) = 0. 



Furthermore, the matrix G* {t, x; y) satisfies 



5G*(t,x;y) ^ d (^*(^^^.y)^^(^)) ^i^^(^) ^ o_ 



dt 



dxi 



Multiplying (2.11) by Tq on the left and on the right, and using the commutation relations 

roAs = Asro, roD^ = -D^ro. (2.12) 



we deduce ( 2.10 ). Then, since Tq = /, (2^) may be re- written as 

vf'^($, xo) = I G(r, xo + z)x(z)e'q--/^x(z )e-^q"'/^G'^ (T, xq + ey; zOro/(q)So (y) "^^"^^^^/"^ 

dydq 



Qf (T, xo + q)Qf (T, xq + ey; q)/(q)roSo(y) 
and (12:81) follows. □ 



(2.13) 



The following lemma allows us to drop the term of order e in the argument of in (|2.8| ). 

Lemma 2.3 Let us define vP(^;xo) = / e^'^'^^-y^ (T, xq, k)roSo(y)^^. There exists a de- 

J (2vr)'* 

terministic function C{e,6) so that 



sup||vf'«-vf'«||^. <C(e,J) 



(2.14) 



and C{e, 5) ^ as /C^ 9 {e, 6) 0. 



The proof of Lemma |2.3| is presented in Appendix^. Note that Theorem 1.1 may be applied directly 
to the moment 

vf(^;xo) = I e^'^-^Ty/(r,xo,k)roSo(k)^ 



and the conclusion of Theorem 2.1 follows. 
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3 The high frequency analysis 



In this section we study the deterministic high-frequency behavior of the Wigner transform and 
estimate the error between the Wigner transform and the solution of the Liouville equations. It is 



weU known |14, 17, ^] that in the high frequency regime the weak limit of the Wigner transform as 
e — > and 5 > is fixed, is described by the classical Liouville equations in the phase space. Here, we 
do not pass to the limit e — > at 5 fixed but rather control the error introduced by the semi-classical 
approximation. As explained in the introduction, this is possible because we are dealing with the 



Wigner transform of a mixture of states that may have strong limits [17] rather than the Wigner 
transform of pure states, which converges only weakly. 

3.1 Convergence on the initial data 



We first show that the assumptions on the convergence of the initial data in Theorem 1 . 1 are purely 
academic, and in particular are satisfied in the time-reversal application. We note that the L^-norm 
of a pure Wigner transform We{t, x, k; q) of a single wave function, such as ( p. 61) , blows up as e — > 
in L^(M'^), because 

\\Wsit;Oh^R^^R^) := I Tr[I^,(i,x,k;C)i^*(t,x,k;C)](ixdk= (2^6)-'^/2||Q-(t,q)||2,^^,^ 



Therefore ( |1.15| ) may not hold for a pure state. Two examples when assumption (1.15) holds are 



given by the following lemma, which may be verified by a straightforward calculation. The first one 
arises when the initial data is random, and the second comes from the time-reversal application. 

Lemma 3.1 Assumption is satisfied in the following two cases: 

(1) Statistical averaging: the initial data is Vq(x;C) = V'(x)y(x/e; C), where V^(y;C) is a mean zero, 
scalar spatially homogeneous random process with a rapidly decaying two-point correlation function 
R{z): E{V{y)V{y + z)} = f ViyXWiy + z;Odfi{C) = R{z) e L^{R% and G C7,(R'^). The 
limit Wigner distribution is given by Wo(x, k) = \ip{x)\'^R(k), where R(k) is the inverse Fourier 
transform of R{y). 

(2) Smoothing of oscillations: the initial data is Vq(x;(^) = '(/'(x)e*^'-^/'^, where ^(x) E Cc(M'^). 
The measure fj, is d^(C) = diOd-d C ^ ^'^j '^^^ 9 ^ L^(]R'^). The limit Wigner distribution is 

T^0(x,k) = |V(x)|2^^(k). 

Proof. We only verify case (2), the other case being similar: 



W^°(x,k)= J e^k.y^(x-^)^*(x+^)5(y) 



so that 



(27r)" 

\\w^-Wo\\l = I (^(x-f)^*(x+f)-|^(x)p)^^(y)l^|^ = fieiymy)?^,. (3.2) 

However, we have /[/^(y)! < 4||'(/'||^4 and 



4(y) = / (^(x-^)^*(x + ^)-|v^(x)|2)'dx^o 

as e — > since ■0 G Cc(K'^), pointwise in y. Therefore \\W^ — Wq\\2 — > by the Lebesgue dominated 
convergence theorem. □ 

Note that if g{C,) and 'ip in part (2) of Lemma |3.1| are sufficiently regular, then ||Wq — VFolb = 0{e) 
so that one may get the order of convergence in ( |1.15D . 
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3.2 Approximation by the Liouville equations 



We now estimate directly the error between the mixed Wigner transform and its semi-classical 
approximation. The dispersion matrix Pq(x, k) = ics{'x.)kjD^ may be diagonalized as 



2 2 

iPo'(x,k) = j;A^(x,k)n,(k), ^n,(k) =/. 

q=0 q=0 



(3.3) 



Here Ilq is the projection matrix onto the eigenspace corresponding to the eigenvalue Ag. Notice that 



the eigenspaces are independent of the spatial position x, hence of the parameter 6; see ( |1.7D -(p^ 

As we have mentioned before, for a fixed 5 > the Wigner transform W^{t,x,'k) converges 
weakly as e — > to its semi-classical limit [/"^(tjX, k) given by 



C/^(t,x,k) = j;n^(t,x,k)n,(k). 



(3.4) 



The functions satisfy the Liouville equations 



'dt 



(3.5) 



with initial data iig(0,x, k) = TrnqWo(x, k)]!^. Our goal is to estimate the difference between 



and in L^{1 
where 



. Let us denote by 7^(x, k) the largest eigenvalue of the matrix (F^F^*)^/^, 



,2\<5 



dkidxj 



H 

\ dxidxj 



dkidkj 
dxidkj / 



Note that ( 1.12 ) implies that 7^(x, k) = 72(x, k) = 7 (x, k), while 70 = 0. The initial data is 
supported on a compact set 5 because Wq is (see (P^)). Then the set 

5 = 



y supp ng(t,x,k) 

t>0, 5g(0,1] 



is bounded because the speed cs{x) is uniformly bounded from above and below for 6 sufficiently 
small ( [1.121) . Therefore we have 7'^(x, k) < 0/5^^"^ with a deterministic constant C > 0. We denote 
75 = sup^ 7''(x, k). We have the following approximation theorem. 

Theorem 3.2 Let the acoustic speed C5(x) be of the form (LS) and satisfy assumptions ( 1.1^ ). We 
assume that the Wigner transform satisfies ^1.15 ) and that holds. Moreover, 

that the initial limit Wigner transform Wq is of the form 



we assume 



T^o(x,k) = ^^x°(x,k)n,(k). 



(3.6) 



Let [/ (i,x, k) = Up(t, X, k)np(k), where the functions Up satisfy the Liouville equations (3.5) 
with initial data u|j(x, k). Then we have 



\\W^it,x,k)-U\t,xM\2<Ci5)[e\\Wo\\me'^'' + e'\\Wo\\m^^^^^^ (3-7) 

where C{5) is a rational function of 5 with deterministic coefficients that may depend on the constant 
C > in the bound ( l.ldj) on the support of Wq . 
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Theorem 3^ shows that the semi-classical approximation is valid for times T <C llnel/75. This 
is reminiscent of the Ehrenfest time of validity of the semi-classical approximation in quantum 
mechanics, see |5|, |9| for recent mathematical results in this direction for the Schrodinger operators. 
The pre-factor constants on the right side of (3.7) are not optimal but sufficient for the purposes of 
our analysis. 

The assumption that initially Wq has no terms of the form npllq with p ^ q \s necessary in 
general for the Liouville equation to provide an approximation to in the strong sense. This may 
be seen on the simple example of the solution 



of the wave equation 



utt 



C Uxx 



with a constant speed c. The cross-terms in the Wigner distribution We = [|ap + \b\^ + ah* e'"^^"^/^ + 
a*6e^*'^*/^]5(k — q) vanish only in the weak sense as a function of t but not strongly. 

The Wigner distribution that arises in the time-reversal application has an initial data that is 
described by part (2) of Lemma |3.1| : 

T^o(x,k) = |x(x)p/(k)/ 



and satisfies the assumption ( |3.6| ) with iig(x, k) 



P/(k) for all eigenspaces because of the 
second equation in (3^). The error introduced by the replacement of the initial data in ( p.7| ) in that 



lx(x) 



case is given by ( |3.2[ ) and is 0{e) provided that x and / are sufficiently regular. 

The proof of Theorem bI3 is quite straightforward though tedious. We first obtain the evolution 



equation for in Section 3.3, and show that it preserves the L -norm. This allows us to replace 
the initial data in the equation for by Wq at the expense of the last term in (3/?). We obtain the 



Liouville equations (|3.5|) in Section pA 



and estimate the right side of (pj) in terms of the iif -norm 



of its solution. Finally, in Appendix |C| we obtain the necessary estimates for the solution of the 
Liouville equation. 



3.3 The evolution equation for the Wigner transform 

The L^-norm of the Wigner transform W{t, x, k; Q) of a pure state, or a fixed Q, is preserved in time 



as follows from the preservation of the L^-norm of solutions of ( |1.5D . We obtain now an evolution 
equation for the Wigner transform of mixed states and show that its L^-norm is also preserved. 
It is convenient to define the skew-symmetric matrix symbol 



P,^(x,k)=Po^(x,k)+ePf(x), 



where Pq is defined by (L7) and the symbol depends only on x: 



Id As, 



Pf(x)=A,(x)Z)^-^(x)-- 



1 dcs 



The latter equality follows from ( |1.6| ) and calculations of the form 

^ 9^-(") = 5^^")"^■®"^+^• 



(3.8) 



(3.9) 



(3.10) 



The following lemma describes the evolution of the Wigner transform W^- 
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Lemma 3.3 The Wigner transform VF|f(t,x, k) satisfies the evolution equation 



'~di 



' + dw^ = 



with initial data W^(0,x, k). The operator C,\ is given by 



4/(x,k) 



P^y, q)e*'^/(z, P) - /(z, p)e-^<^P,^(y, q 



dzdpdydq 



(3.11) 



(3.12) 



(7re)2'i ' 

where 0(x, z, k, p, y, q) = 7((p — k) • ^/ + (q — p) • x + (k — q) • z). The integral of the trace and the 
L^-norm of the Wigner transform are preserved: 



and 



J TrWf(t,x,k)dxdk = J TVWf(0,x,k)dxdk 
j Tr[Wf(t,x,k)Wf (t,x,k)]dxdk = j Tr[Wf(0,x,k)Wf (0,x,k)]dxdk. 



(3.13) 



(3.14) 



This lemma is verified by a direct calculation that we present for the convenience of the reader in 
Appendix M. 



Note that the solution of (3.11) with self-adjoint initial data remains self-adjoint and the L^-norm 
is preserved. Therefore, we have the following corollary. 



Corollary 3.4 Let lVo(x, k) be a strong limit o/ We(0) in L?' , which exists by assumption 

Then the solutions W^{t, x, k) and W^{t, x, k) of \3. 1 ^ ) with initial conditions VF|P(x, k) and Wo(x, k), 

respectively, satisfy 



\w!{t)-wlyzn\2 



\W^-Wq\\2^Q ase^Q. 



This shows that in the analysis of ( |3.11| ), we can replace strongly converging initial conditions by 
their limit, and consider then the limit of W£(t,x, k) as e — > with fixed initial conditions. This is 
done in the following section. 



3.4 Derivation of the Liouville equations 

We consider in this section the solution W^(t, x, k) of the evolution equation ( 3.11| ) with fixed initial 
data iyo(x, k) and show that it may be approximated by the solution of the Liouville equation. We 
split the operator £^ = Ce'^ + eCi'^ , where 

fizdpdydq 



4V(x,k) = j (p/(y,q)e^<^/(z,p)-/(z,p)e-^<^P/(y,q) 



(vre) 



2d 



i = o,i, 



and the symbols are given by (|l.7|) and (|3.9| ). The operator is given explicitly by 



4'°/(x,k) = / e^e^-P^-y 



c<5(x - ^-)ipjD^ f{y., p) - C5(x + y)/(x, p)ipjD^ 



dpdy 



+£ / e^(k-p) y (^i^ 



d_ 

dxi 



C5(x- y)/(x,p) 



+ 



d 



c<5(x + y)/(x,p) 



{2'kY 
dpdy 



^e,5 + ^e,5- 



We recast the operator £^^^5 as 

£0V(x,k) = C5(x)[ifc,I)V(x,k) -f{^,k)ik,D^] - fel^V(x,k) + f{^,k)kjD^) 
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with the correction term 



7^o^^/(x,k) = - /e^(k-p)-y (c5(x-^)-C5(x) + ^yVc5(x))ip,Z^V(x,P 



(c<5(x+ y) -cs{x) - |y- Vc5(x)^ f{x,p)ipjD^ 



dpdy 
{2-kY' 



Similarly, we have 



£oy(x,k) = ^-D^^ (C5(x)/(x,k)) + (c,(x)/(x,k)) D^+eUfJ 



d 



e d 



02 



2 dxj 



2dxj 



with 



7^02,/(x, k) = ^ y e^^'^-^yy (i?^- A ^|c,(x - ^) - C5(x)} /(x, p 




5xj 



{c5(x + y)-C5(x)}/(x,p) 



The operator is given explicitly by 



4'V(x,k) 



pi(k-p)-y 



(x - ^)/(x, p) - /(x, p)Pf (x + 



dpdy 

{2-kY 



= (x)/(x,k) - /(x,k)Pf (x) +7^l,,/(x,k) 
with the correction TZ], ^ defined by 



7^,\,/(x, k) = / e^C^-P)-^ (Pf (x - ^) - Pf (x) ) /(x, p) - /(x, p) (Pf (x + ) - P^ (x) 



Putting together the above expressions, we obtain the following equation for 



dW^ 1 r - r 



[^r,^p/l 



dt e " " e 
with 7^f = 7^°^5 + nfi^ + 5. Here {/, g} is the standard Poisson bracket 

{/,5} = Vk/-Vx5- Vx/-Vk(7 

and [A, B] = AB — BA is the commutator. We now introduce the expansion 

= U^ + eUi + Ul^. 

We insert this ansatz into ( |3.17 ) and equating like powers of e obtain at the order 

[Po^?7^]=0, 

which is equivalent to 



1 ' 



q=0 



q=0 



(3.15) 



dpdy 

{2ttY' 
(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 
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where = HgU^Ug, and for q = 1,2 one has = uplg with Ug = Trllg. The matrices Ilg are 
projections on the eigenspaces of Pq, as in ( |3.3D . This means that the matrix does not have 
off-diagonal contributions in the eigenbasis of Pq. The equation of order 0{e^) is given by 

-Qf = [Uf,P^] + [U',Pf] + Yi(iU',P^} - {PIU'}). (3.21) 
Multiplying the above equation on both sides by Ilg yields 
r)T\ TI^T\ 1 

-S^=Ug[U\Pf]ng+-Ug{{U^,P^}-{P^,U'})Ug. (3.22) 



This is nothing but equation (6.16) of reference |14] for the Wigner matrix without consideration of 
mixtures of states. The only difference is that the leading order term Pq depends on the parameter 
6. This, of course, does not change the algebra, and following [14] one obtains a system of decoupled 

dui 



Liouville equations for = TrC/^ , g = 1, 2, 



dt 



1=0 (3.23) 



with initial data ^^(0, x, k) = Tr[n<y(k)Wo(x, k)ng(k)]. The zero eigenvalue component of the matrix 
, that is, C/o(t,x, k) = no(k)Wo(x, k)no(k), does not change in time. 

We have to show that the terms and in ( |3.18| ) are small. In order to uniquely characterize 
f/f , we assume that it is orthogonal to the terms of the form ( 3.2C| ), that is, 

c/f = ^np[/fn,. (3.24) 

Then, (|]20|) and imply that 

li^Uilip = ^ IlmB{U^)Iip, (3.25) 

where 

B{U') = [U\Pt] + ^({^^^'o'} - {pLu'}). 

We now analyze the term C/|g in ( |3.1g| ) and show that it vanishes in the limit e ^ 0. The equation 
for the C/l e is 

-^ = -Ciul + Se, (3.26) 

where 

r<5 D^i I ^ (<TT& d5t rD<5 TT5^\\ ^^^1 T>StTT5 I ^tt5\ 



Se = e[[UlPt] + 1^{{UIP^} - {PL Ul})) -^-Qf- KiU' + eU^). (3.27) 

The initial condition for (|3.26|) is C/|£(0,x, k) = — e[/f(0,x, k) because of (|3.6D , which implies that 
Wo(0,x,k) = C/'5(0,x,k). We now use the fact that is skew-symmetric to obtain the bound 

WlMh < 4Uf{0)\\2 + f \\Se{s)hds. (3.28) 

The analysis of the convergence of the difference of and to zero thus relies on estimating 
the error term S^- The relevant bounds are provided by the following two lemmas. Here we denote 
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Lemma 3.5 There exists a constant C > that depends on the constant in the bound ( 1.16 ) on the 
support of Wo, and on the constants Di in ( 1.1^ so that 



Seh < C 



(3.29) 



Lemma 3.6 The H'( 



-norm, s = 1,2,3, ofU^{t) is bounded by 



(3.30) 



Here = TrpqWoHq], the initial data for the Liouville equation (3.2c), the constant Cg is a 
deterministic rational function of 6. 



Note that the prefactors of the type (5"™ in Lemma 3.5 are not as important as the terms 
since the latter grow exponentially in js ~ d~^^'^ according to Lemma 3.6 . 
Proof of Lemma |3.5| . Observe that thanks to ( p. 27 ), we have 



< c 



\Uf\\m+e 



dt 



+ \\ni{u^ + £ui) 



(3.31) 



We have the following bound for 



<5. 



1 WH- 



(3.32) 



with a constant C > that depends only on the constant in the bound ( |1.16| ) on the support of Wq 
and on the constants Di in ( 1.12| ). Indeed, expression ( 3.24 ) implies that ||C/f < C52~'^\\B{U^)\\h^ , 
while we have \\B{U^)\\h^ < C6-'-2\\U^\\ fjs+i so that (|3.32| ) follows. This bound is by no means 
optimal but will be sufficient for our purposes. Furthermore, we have 



dt 



< C 



B 



In order to complete the bound ( |3.29| ) for we show that 



Ce 



(3.33) 



(3.34) 



We only consider TZ^^^ as the corresponding bounds for the operators 'R-^'^g and TZ^ g are obtained 
similarly. We split 7^°^5 as TZ^^^g = Iqi - Hoi- We have 

1 



loif 



1 



^(i'-P)-y ( cs{^ - ^) - q(x) + |y • Vc5(x)) ip,D^f{^, p)^ 



4^/ 



d'^csix 



dxidxm 



^ 'iPjD^fi^^P)7^^ds = ^l {£ - s)Ioi{s)fds. 



(27r)« 



4e 
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Moreover, we obtain that 



dxidxm. 



dxvdxm' 



xpjqrD^f{x,p)r{^,q)D'' 



Tr / e'^'i ^^'^yiymyi'ym' 



. dpdydcidzdxdk 

a2c5(x-^)52c5(x 



2 ' 



dxidxr, 



dxi'dXm' 



.p,(?,Z^V(x,p)r(x,q)D^ 



. dpdydqdx 
(2^ 



C 



Therefore the Minkowski inequahty imphes that ||/oi/||2 < Ce5 '^^^ ||%/||//2, and the same 
bound holds for //qi. The operators TZ^'^ and TZ] may be bounded in a similar way as ||7^e^/||L2 + 
WR-lfh^ < CeS^^/^\\f\\H2. Therefore we have the bound ( Pl and then (|09|) follows from (|Ol| )- 

dJl). □ 

Theorem 3^ now follows from the bound ( 3.32| ) for Uf, the bound ( 3.28| ) for f7|g, and Lemmas 
|3.6| and 3.5. It only remains to prove Lemma which is done in Appendix O. 



4 The Liouville equations in a random medium 

We formulate in this section the main result concerning the convergence of the expectation of the 
solution of the Liouville equation (|l.9| ) to the solution of the phase space diffusion equation ( |1.10 ) 
in the limit 5^0. We also show that the values of the solution of the Liouville equation at 
different points in the phase space become independent in this limit. This allows us to establish the 
self- aver aging property in Theorem l.L 



4.1 Preliminaries 

We let Cm ■■= C([0,+oo); (M'^)™), and for any iii, • • • , > we denote by Cm{Ri,--- ,Rm) ■= 
C (^[0, +oo);5^~^ X ••• X Sf^^, where 5^"^ is the sphere in W'^ of radius R > centered at 0. 

We also let vr^ : Cm. (R'^)"^, t > 0, he the canonical mapping TTt{K) = {Ki{t),--- ,Km{t)), 
K = (Ki,--- ,Km) G Cm. For any u < f we denote by M.^ be the cr-algebra of subsets of 
Cm generated by vr^, t G [tt, f], and let M.m '■= and Tm be the filtered measurable space 

Cm,Mm, (■^"'*)t>o)' ^'^^ ^ ^ ^^^''^ denote C{A) := cj{ci(x) : x G A}. 

We suppose further that $7 — > M is a scalar , measurable, strictly stationary, zero mean 

random field that satisfies assumptions presented in Section L2, that is, it satisfies the almost sure 
bounds ( 1.12| ), is exponentially (/)-mixing ( 1.13| ), and has a C°°-correlation function /2(x). 
We define the differential operator 

d d 

£F(k) = |k|2l),,,(k)a|^ ,^F(k) + |k|ii;p(k)9,,i^(k), F e Co°°(M'^ \ {0}) (4.1) 

p,q=l p=l 

with the diffusion matrix D given by ( |l.ll| ) and the drift E defined by 

d 

^p(k) = -coE / «5xp,x„x, ^(cosk)ds, Vp = l,---,ci. 
9=1 n 
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A simple calculation shows that £ is a generator of a diffusion on 5^^ ^ given by Ito S.D.E. 



I dk{t) = \k{t)\ (E{k{t)) dt + V2DV2(fc(t)) dB{t) 
\ k{0) = ko / 0. 

Here E = (Ei, • • • , E^i) and B(-) is a d-dimensional standard Brownian motion. 



(4.2) 



Remark 4.1 A simple calculation shows that the diffusion fc(-) given by (|4.2| ) is symmetric. Indeed 
the generator can be written in the form 

d 

£F{k) = dk, {\k\'Dp,,{k)dk,F{k)) , F e C^iR'' \ {0}). 

p,q=l 

For any k / we denote by Ok the law of such a diffusion starting at k, which is supported in 
Ci{k), k = |k|. 

Remark 4.2 The matrix D := [Dp^g] is degenerate since D(k)k = for all k € M'^ \ {0}. It can be 
shown however that under fairly general assumptions its rank equals d—1. 

Proposition 4.3 Suppose that R{0) > 0. Then, the rank o/D equals d—1. 

Proof. Suppose that cq = 1 and let i^k := [p G M*^ : p • k = 0] be the hyperplane orthogonal to 
k. Then, 

D^iik) = -1 Pjl^^,^R{sk)ds = (^l^^e''^-PpmPiR{p)dp^ ds 

" ¥^ j^JrnPlR{p)dp 

and hence for any ^ G M*^ we have 

(D(k)^, ^) = Dmi{k)U^i = (p • ifR{v)dp. (4.3) 



Suppose that ^ G H\f^. Then, since i?(p) > the left hand side of (O) is nonnegative. We 
claim that in fact (D(k)^,^) > 0. Indeed, if otherwise then, since R is continuous, we would have 
i?(p)(p • ^)2 = for all p G i^kj which is impossible due to the fact that -R(O) > and the set 
n -ffk has the linear dimension d — 2. □ 

The above argument shows that D(k) is of rank d— 1 if there exists po G i/k such that -R(po) > 0. 
On the other hand, if -R(p) = for all p in the plane i^k then D(k)^ = for all ^ G W^. Therefore 
the matrix D(k) either has rank d—1, or vanishes identically. Another condition ensuring the latter 
does not happen is the radial symmetry of R{-). 

4.2 Two particle model 

We would like to show that solution u^{t, x, k) of (|1.9| ) decorrelates in the limit 5 ^ at two different 
points, that is, that 

E|u'^(t,xi,ki)n^(t,X2,k2)} -E|n^(t,xi,ki)}E{n^(t,X2,k2)} ^0 as 5 ^ (4.4) 
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provided that ki 7^ k2. Recall that u^{t,x,'k) may be represented as 



4(.T, X, k) = <(x^(r, X, k), -K'{T, X, k)), 



where is the initial data for (|1.9D, and 



dX^(t) 
Jt 

di 



^(X^(t),K^(t)), X^(0)=x 

--^(X^(t),K'5(t)), K^(0) = 
ax 



(4.5) 



In order to establish (4.4) we have to consider motion of two particles that may start at the same 
physical point but are moving in different directions. The equations of motion for two particles are 
governed by the Hamiltonian system 



dt 



flf — V x^g I •^m V'') ■'•-mi '^mj: '^m ■^mi 



S^m'' (0; Xm, k^) — Xmj f^rri {0] , lif^ 



We will assume that xi = X2 = 0, and 

ki / 0, ka / and ki / kg- 
The above system can be rewritten in the form 



m 

m = 1,2 



(4.6) 



(4.7) 



dt 



di 



Co 



k'-^^f-x k 



(4.8) 



l^m (i)Xm,km)| 

^ i'^m^ (0; Xjji, km) = 0, fcm'' (0; x^, k^) = k^,, th = 1,2. 
The main result of this section is the following. 

Theorem 4.4 Suppose that the random field ci(-) satisfies the assumptions in Section \L.^ and that 
d > 3. Then, the laws of processes 

{k[^\-),x[^\-),k^2^\-)^X2\-)) determined by (pj), converge 
weakly inC^, as 5 — > 0, to the law of {ki{-),xi{-),k2{-),X2{-)), where kj{-), j = 1,2 are independent 
symmetric diffusions given by (4^) starting ai kj, j = 1,2 respectively and 



Xjit) 



-Co 



z 

j kj{s)ds, 



J = l,2. 



Theorem LI is a simple corollary of Theorems 3^ and 4^ 
Proof of Theorem ll.ll. First we observe that 



2 

j j iW^it,x,k) -U\t,x,k)S{k)d-k dx<\\S\\l2 j |W/(t,x,k) - f/^(t,x,k)|2dk(ix ^ 

as {e,5) — > in /C^ and this convergence is uniform in realizations of the random medium provided 
that the bounds ( [1.12 ) are satisfied. Therefore it suffices to study s^{x) = J C/^(t, x, k)S'(k)(ik. We 
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observe that 



E 



is'^(x) -s(x)fdx 




{U^t, X, k) - W{t, X, k))5(k)dk 



dx 



E<{ 5*(k)i / (C/'^*(t,X,ki) - l^*(t,X,ki))(C/''(t,X,k2) - VF(t,X,k2))5(k2)dkidk2dx 



with s(x) and VF(t,x, k) as in the formulation of Theorem lA. Theorem iA imphes that 

E{c/'^(t,x,k)} ^ T^(t,x,k), E{f/'5(t,x,ki)[/'5(t,x,k2)} ^ VF(t,x,ki)T4^(t,x,k2) 

pointwise in x and k. Recall that the functions ?7'^(t,x, k) and x, k) are uniformly compactly 
supported and bounded in L°° . Therefore the Lebesgue dominated convergence implies that 



E 



II ||5^(x)-.-(x)fdx|^0 



and the proof of Theorem |l.l| is complete. □ 



5 Proof of Theorem 14.4 



Before we present the proof of this result we wish to spend a few words to lay out its main ideas. 
They are based in large part on the ideas of |15] where the phase space diffusion equation for the 
limit of the expectation of the solution of the Liouville equation with the Hamiltonian ff''(x, k) = 
+ \^V{x/6) has been obtained. The two-particle case introduces some additional difficulties 
into the problem. Our first step in the proof, in Section 5A below, is to replace the processes 
, k2{-)) by (ii(-)j ^2(')) that agree with ^2(')) up to certain stopping times. These times 

are determined by the stopping rules, introduced by multiplying the Hamiltonian A'^(x, k) by several 
cut-off functions. Their role is to prevent the trajectory of each particle to self-intersect and also 
not to allow the particles to get too close to each other. We shall prove tightness of such modified 
processes by showing that for any bounded, positive and continuous function F one can find a 
constant C > such that F{lf{t),l2{t)) + Ct, t > are sub-martingales (see e.g. 



Theorem 



1.4.6), cf ( ^.291) . This fact will be established thanks to the decorrelation properties of the random 



field VxCi(-). More precisely, the latter imply mixing lemmas contained in Section 5.2. The second 
ingredient of the proof is a perturbative argument that allows us to replace the trajectory x^f\-) (in 

fact its modification y[^\-) that arises from the replacement of by l^) by a linear approximation 
over the time interval that is much longer than the correlation time (that we recall is of order 0(6)) 
yet is sufficiently short so we can control the accuracy of the approximation, cf. Lemma ^.4| . In 
order to ensure that the approximate motion (under linear approximation) is not transverse to the 
direction of the field at a given time, which could prevent us from using the decorrelation properties 
of the field, but is rather propelled forward, we have to introduce another stopping time rule, cf. the 
condition on the scalar product of wave number directions contained in (|5.5| ). 

Conducting the proof of tightness we also identify a certain martingale property of any limiting 
law of (ii(-)) '2('))i as 5 ^ that holds up to the aforementioned stopping time. By proving that this 
time goes to infinity with the removal of the cut-offs we are able to prove both the weak convergence 
of the laws of {kf(-), fc2(')) s-'^d identify a well-posed martingale problem associated with the limiting 
measure. This step is done in Section 

With no loss of generality we shall assume throughout this section that cq = 1. 
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5.1 The cut-off functions 

Let p,q > and > be integers. Let M be chosen in such a way that 

M>|ki|V|k2| and |ki| A |k2| > 



Let ki 7^ k2 be such as in the statement of Theorem 4.4. Denote 

1 



K^:= (k,k'): (k,ki)jj, >1 



, (k',k2)K. > 1 



1 



N + 1' 

and choose a positive integer such that 

7JV := inf ||k-k'| : (k,k') G Ktv} > 0, 



iV + 1 



(5.1) 
(5.2) 

(5.3) 



that is, the cones of aperture l/(iV + 1) centered at ki and k2 are separated. As a consequence of 
(|5.3|) we may choose a positive integer q so that 



A7v(p) := inf 



k-pk' 



A 



-k' - pk 

P 



0, 



1 



, (k,k') GKjvl> > -. (5.4) 
Q 



We define now several auxihary functions that will be used to introduce the cut-offs in the dynamics. 
The function i/j-.R'^ x {Sf'^)"^ [0, 1] is C°° and has the property that 



^(k,li,l2 



if k • li > 1 

and 
if ic • li < 1 



1 



and k • I9 > 1 



N+l 

< |k| < M 

2 



N+l 



N+l 

or |k| < (2A/)-i 



or k • I9 < 1 



or 



k| > 2M. 



(5.5) 



The function 



X Ci ^ [0, 1] is for a fixed path K{t) and satisfies 



(/>fc(y;K) = < 



1, if inf 

o<t<t[-\ 

0, if inf 



y-J Kis)ds 


y-jK{s)ds 



> 



< i. 



Here t^^' := kp~^ and, by convention, K{s) := K{0), s < 0. The function : M*^ x M"^ x C2 
is smooth when the paths Ki{-), i^2(") G Ci are fixed. We let 

pi :=2«[8(l + Z)o)b 
and s^^^ := kp^^ be a sub-partition of t^, and define 



(5.6) 

[0,1] 
(5.7) 



6(yi,y2;i^i(-),^2(-)) = { 



if inf 



o<t< 



,{pi) 



and 



inf 



0<i< 



,{pi) 



0, 



if inf 



0<i< 



,(pi) 



or 



inf 



0<t<s 



(Pi) 



yi - / K2{s)ds 


t 

y2 - / Ki{s)ds 


t 

yi - / K2{s)ds 


y2 - / Ki{s)ds 




> 



> 



< 



< 



(5.8) 
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For j = 1,2 we set 



(p) 



(5.9) 



1, if < t < 

^ </>fc(y;K(.)), if4^)<t<tg,. 

Each $j(-) shall be used to modify the dynamics of the corresponding particle in order to avoid a 
possibility of self- intersections of its trajectory. The cut-off function 

for t G [4^\4+i) and A: > 1 

'(/'(k,i^(0),K(0)) for t G [0,tS^^) 

will allow us to control the direction of the particle motion over each interval of the partition as 
well as not to allow the trajectory to escape to the regions where the change of velocity can be 
uncontrollable. The cut-off 



(5.10) 



H(t,yi,y2;/^i(-),^2(-)) 



1, 



if < t < t 

,(pi) 



(p) 
1 



^ 6(yi,y2;/^i(-),^2(-)), if < t < 4+1 ^nd < 4 

is introduced in order not to allow the two trajectories to come too close to each other. Note that 
this cut-off is "switched on" only after time t = t\ so as to allow the two particles to separate 
initially. After this time it is updated every 1/pi time step, that is, more frequently that the cut-offs 
that control the self-intersections of each trajectory that are updated only at each \/p time step. 

The following lemma can be checked by a direct calculation. Both here and in what follows we 
denote by D,^p the partial with respect to the (3 component of the given vector variable. 

Lemma 5.1 Let m = (mi, • • • ,771^) he a multi-index with nonnegative integer valued components, 
d 

m = rup. There exist constants C3, C4 > depending only on AI,N,p,q, m such that 



.(pi) 



(5.11) 



\D^^j{t,y)\ < C3, \D^B{t,yi,y2)\ < C4, j = 1,2. 
Let K = {Ki,K2) G C2 and denote 

ej{s,yi,y2,l;K) := ^{s,l;Kj)^j {s,yj; Kj) E {s,yi,y2; K) , 



(5.12) 



^j(.s, yi,y2,y'i,y2,hK) := @j{s, yi,y2, 1; K)ej{s,y[, y2,l; K). (5.13) 
We also introduce a random transformation of paths A'(-) = (i^i(-), K2{-)) for any K £ C2 given by 



1 + VSci 



K,{t) 



Kj{t), t > 0. 



Finally, let us set 

F,it, yi, y2, 1; K) = e,(t, 6yi,6y2, 1; K)Vy^ci (y,) |1|, j = 1, 2. 
The modified two particle system with the cut-offs that we will consider is given by 



IF 



( dyf\t) 



dlfjt) 

dt 



1 + V^^ 



Cl 



Ij (^i 5 kj) 



{ y;"^(0)=0, if)(0)=k„ j = l,2. 



where the path I 



(l^i\-)-,l2' {')) is obtained from /(•) by the transformation ( ^.14| ). We will 



(5.14) 
(5.15) 



(5.16) 



denote by Q\--M,N,p,q) the law of {if {■),y'{' {■),y)^' {■)) for a given 5 > over C4 
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5.2 The Mixing Lemmas 

For any t > we denote by Tt the fj-algebra generated by (Z^'^^(s),4^^('S))) s < t. Throughout 
this section we assume that Xi,X2 : M x M*^ x M^^ — > M are certain continuous functions, Z is a 
random variable and ^1,52 are M*^- valued random vectors. We suppose further that Z,gi,g2, are 
.T^t-measurable, while Xi,X2 are random fields of the form Xj(x) = Xi (ci(x), VxCi(x), V^ci(x)) , 
satisfy lim ||Xj(x) - Xi(0)||oo = 0, i = l,2. We also let 

|x|^0 

u{ei,e2) ■.= E[Xi{ei)X2{92)], {61,62) e {R'^f. 

The following mixing lemmas will be of crucial importance for us in the sequel. 
Lemma 5.2 Assume that r,t>0 and 



(5.17) 



inf 

u<t 



yf{u) 



r 

> - 
- 6 



-a.s. on the set Z for i,j = 1, 2. Then, we have 

\E[Xi{gi)X2{g2)Z] -E[U{gi,g2)Z]\ < 2(P (^^) ||Xi|U||X2||oo||^||i. 



(5.18) 



(5.19) 



Proof. The proof is a modification of the proof of Lemma 2 of |15] so we only highlight its main 
points. Choose an arbitrary rj > 0. By a suitable modification of gi,g2 on the event Z = 0, so that 
the modified r.v. remain .Fj-measurable, we can guarantee that ( 5.18 ) holds P-a.s. Let 

i = (zi, • • • ,i,) G and C; := [h/2^\{h + l)/2^'n x • • • x [id/2^'\{id + l)/2''^) 

and 

Q := ((2ii + l)/2*-^i+\ • • • , {2i, + l)/2*^i+i). 
Here Mi > is a sufficiently large integer so that 

||Xi(x)-Xi(ci)||oo <??, ViGZ^xGCi, i = l,2 (5.20) 

and 2"^^i < r/(20(5). We let 

Aj := [z : dist{z, d U Cj) > r(2(5)-^] 

and 



2) 



{yf\s),yr{s)):s<t 



Let us denote by li j the indicator of the event [{gi ,52) £ CixCj] and the event yli j = [to : 2)|^^ (lo) C 
A,j]- Note that 



E [Xi (50X2(52)^] = J^E [Xi(5i)X2(52)^/ijXA,J . 



(5.21) 



Using precisely the same argument as in 1 15 1 we prove that ZI^ jXAi^ is C(I?ij)-measurable for each 
i,j G 7j'^. Note however that the right hand side of ( 5.21 ) is equal, up to a term of order 0{rj), to 

X;E[Xi(ci)X2(cj)Z/i,jXA,J. (5.22) 
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The random variable Xi(ci)X2{cj) is however C(Ci U Cj)-measurable. Therefore we can write, see 
e.g. |] p. 171, that 



|E [Xi(ci)X2(cj)Z/ijXA,j] - t/(ci,Cj)E [Z/ijXA,j] 



(5.23) 



1.J 



< 



1.J 



|E[Z/ijXA,j]| ll^l||oo||^2||c 



However, U (q, cj) equals, up to a term of order 0{r]), to U (51, 32) on the event corresponding to /jj. 
The conclusion of Lemma |5.2| follows upon the passage to the limit Mi — > +00 and t] I 0. □ 



Lemma 5.3 Assume that r,t are as in the previous lemma. Let KXi = 0. Furthermore, we assume 
that §2 satisfies ( 5.18| ), 

{6), N 

r + ri 



inf 

u<t 



91 



yf{u) 



> 



-, i = i,2 



and 



I I ^ ^1 

|9i - 92] > -J, 



(5.24) 
(5.25) 



for some ri > 0, ¥-a.s. on the event Z ^ 0. Then we have 



\E[Xi{gi)X2{g2)Z]-E[U{gi,g2)Z]\ < C^c^^/^ (^) ^^/^ (^) ||Xi||oo||X2||oo||^||i (5.26) 
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for some absolute constant C5 > Here the function U is given 



Proof. We prove that the left hand side of (5.26) is bounded by 

C6(A(^) ||Xi||oo||X2||oo||^||i. 



(5.27) 



This together with the result of the previous lemma imply (5.26). 

Let T] > and Mi be as in the proof of Lemma |5.2| , and in addition 2"*^^ < ri/(205). Note that 
X2(cj)Z/ijXyii j (in the notation of the proof of Lemma |5^ ) is C{Di j U Cj)-measurable. In addition, 
we have dist{Ci, Dij U Cj) > ri{26)~^ thus, using the mixing coefficient as in e.g. [@] p. 171 we can 
estimate 



|E [Xi{ci)X2ici)Zh,jXA,] I < 2cl> (^) ||Xi|U||X2||oo||^||i. 



26 



On the other hand, we have /ij 7^ only if |ci — Cj| > ri{26) ^ , which in turn implies that 



ri 



|C/(Ci,Cj)| <C7(A U7 ll^l||oo||^2||oo, 



with the constant C7 independent of 77 > 0. Summarizing, we have shown that 

ri 



J2\U{cuCi)E [Z/ijXA,j]| <C80(^) ||Xi|U||X2||oo||^||l, 



1,J 



with the constant Cg independent of > 0. Letting i] ^ and using (|]2D we conclude (|5]2|). □ 
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5.3 Tightness and the martingale property of hmiting measures 

In this section we prove tightness of the family Q^{-; M, N,p,q), 6 G (0, 1] and show that any weak 
Umit point Q{-;M,N,p,q) of this family as 5 — > 0, has a certain martingale property. 

Let C^'^^^^P^i be a random partial differential operator defined on C^((M'^)^) as follows. For any 
K = {Ki,K2) G C2 and G G C^{{R'^f) we set Y = (^1,^2) G C2, 



Yi{t) = / Ki{s)ds, 



1,2, 



(5.28) 



where 

<^i4t;Ki) ■.= <^i{t,Yi{t);Ki), ^,{t;Ki) := ^{t, K,{ty, Ki), E,{t;K) := E{t,Y^{t),Y2{t); K) . 
We let 



M,N,p,q 



G)(ki, ks; K) := e?(t; K)£^,G{ki,k2) + OUt; K)£k.G(ki, ks), 



with /^k,, i = 1, 2 given by (4.1). 

Let C £ Cb((M'^)^"') be an arbitrary nonnegative function, let < ti < • • • < t„ < t < u and 
define C{K) := C{K{ti), ■■■ , K{tn)). We will show that for any function G G C^HW^f) there exists 
a deterministic constant Cg > such that 



{[G{i?iu)4'\u))-Gii?{t)4'\t))] ai['\-)4'\-))} 



(5.29) 



<C,{u-t)E[Cil['\.),li'\.))], VC(-),'^G(0,1]. 



The choice of the constant Cg may depend on a particular function G but should be the same for all 
the spatial translates of G, and may not depend on the test function C. This, according to Theorem 
L4.6 of m, implies tightness of the laws of 4i\-) , , 6 G (0, 1] over C2. 

Additionally, we prove that if Q{-; M, N,p, q) is any limiting law of Q^"{-; M, N,p, q), as (5„ — > 
then 



lim E' 



[[G{lf-\u),li'-\u))-G{lf"\t)4'"\t))] C(4'"^(-),4'"^(-))} (5.30) 



M,N,p,q 



G){K{s)-K)ds 



CiK)\QidK;M,N,p,, 



for any u > t. This property will be used in the next section to identify the limiting law of 
(fcf)(-),fcf (•)), as^^O. 

Throughout the remainder of this section we suppress writing both the superscript 5 and the 
cut-off parameters M, N,p, q of the respective measures. With no loss of generality we assume that 
there exists ki such that s^^^^ < t < u < s^^^^i, cf. ( |5.7| ). Given s > o" > 0, we define the linear 



approximation 
and 



Lj{a, s) := yj{a) + (s - (7)ij{a), 



Rj{v,a,s) := (1 - v)Lj{a,s) + vyj{s), j = 1,2. 
The following simple lemma can be verified by a direct calculation. 
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Lemma 5.4 Suppose that s > a. Then, 

\y^{s)-L,{a,s)\<^^^^^^^+DoV6{s-a), V5>0, j = l,2. 

Remark 5.5 Throughout this argument we use 

a{s) := max[t, s - S^''^'] for some 71 £ (0, 1/8). (5.31) 

The above lemma proves that for this choice of a the linear approximation Lj{a,s) of the particle 
position given by yj{s) is exact, up to a term of order 

0((53/2-27i). (5.32) 

We begin now the proof of ( ^.291) . Our strategy is based on the perturbation method: the 
trajectory is approximated by the iterated hnear approximation sufficiently many times so that 
the error becomes deterministicaUy small. The terms that involve the linear approximation are 
potentially large but are handled with the help of the mixing lemmas. Note that 

Gil['\u),li'\u))-G{l['\t),li'\t)) (5.33) 



We can rewrite ( ^.33| ) in the form 
where 



71 E/ A,,aG(Zi(s),Z2(5))F,,a (^S,- 



-,lj{s) ds. 



/(i)+/{2)+j(3)^ (5.34) 



u 

u s 

^^'^ -■= -sT.il [ [ D,^,pD,^,^G{h{p)Mp)) 

X Fj,a I s, — — ) Fi,i3 [p, — —^,h{p) j dsdp. 

5.3.1 Term E[/(i)C]- 

The term /^^^ can be rewritten in the form 

where 

u 
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and 



u 1 



i,a i,0 t 

X (yi,/3(s) - Li^p^a, s)) ds dv. 



(5.35) 



Note that we have replaced by its linearization Lj in the term J^^\ The linear approximation is 
always "propelled forward", which allows us to use Lemma to handle the term E[j(^)C]. Suppose 
that k is such that s, t G i^+i)' I'^call also that s,t,u £ [s^^\ s^^^i), and let us fix one trajectory 
by setting, for instance, j = 1. We will use Lemma 5.2 with Xi{x) = — VxCi(x), X2{x) = 1, 



and gi = Li{a, s)6 ^, cf. ( ^.121) . We need to verify ( 5.18| ). Suppose therefore that Z 0. For p G 
[0,4^ij we have \Li{a,s) - yi{p)\ > (2g)-^ provided that <5 < (2g)"^/(^"'^i). For p G [4^2i,cr] 
we have 



{L^{a,s)-y,{p)).h (41 

> (. - a)h (-) • h (4^2i) + / [l + 

p 



(5.36) 



/i(pi)-/i(4-i)'^Pi 



iv + ly ' 

provided that 6 < 1/Dq. We see from (5.36) that ( |5.18| ) is satisfied with r = ^1 — tyTt) ~ ^'^'^ 
j = l- 

We verify next that gi is also separated from y2{p)S^^, P G [0, o"]. Consider two cases. First, 
when s,t £ [0,4^^), 

using condition ( |5.3| ) we obtain then that there exists 7^ > depending only 

on N such that 
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Z/2(P) 



6 



> 



Suppose then that s,t > 1/p and s,t G [■^i^^^ ■^fcj+i)- Then we have for p G [0,s^^^'], with pi 

given by (5/7), |-Li((T, s) — y2{p)\ > (2(7)"^, provided that 5 is as above. For p G we get, 

thanks to (5.7), 



\Li{a, s) - y2(p)| > Li{a, s) - y^ ( sf^^ 



y2{p) 



> > — > (s — CJ) 1 



iV + 1 



provided that 6 < {Aq)-^^-^^K 
Using Lemma |5.2| we estimate 



< 



MDq 
V6 



||VG||i^((iR.).)E[C] / 0(Cio— 



S — CJ 



ds 



(5.37) 



<Cii(5)(^/-t)||VG||ioo((K.)2)E[C], 
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where Cio := min[7^, 1/2(1 — 2/[N + 1))], and Cii(5) depends only on 5 and vanishes as 5 
On the other hand, the term J^^^ defined by ( |5.35 ) may be written as 



0. 



J(2) = j(2) + j(2)^ 



where 



J- 



(2) ._ 



u 



and 



J, 



(2) ._ 



A5/2 EE 



Mil 



, Ri{9v,a,s) R2{ev,a,s) 
^yun^yuP^j,a s, , ,lj{a} 



§5/2 J J J "yfe.7-yi.M- J," y-' § ' § 



(5.38) 



X Aj,QG(Zi(fT), 12{(t))v ivi^pis) - Li,p{(T, s)){yk,^{s) - Lk,^{a, s)) ds dv dO. 
The second term may be handled easily with the help of Lemma |5.4| and ( |5.32| ). We have 
|E[jf C]| < Ci2M[C]||VG||i^((K.)2)(n- 0^5/253-471 < Ci35i/2-^^H^^-i)E[C]||VG|| 



(2) 

In order to estimate we split it as 



t(2) _ t(2) , x(2) 
•^1 — "^l,! "T "^1,2 



(5.39) 
(5.40) 



where 



J 



(2) ._ 
1,1 ■" 



^EE / / A„.G(Z,(.),/,(a))Z^,,,,F,. (.,^i^,^^,^.(a)) 



d 



{s-Pi)-;—k.i3{pi)dsdpi, 
dpi 



and 



7(2) .. 
"^1,2 •" 



with 



X Ci 



d 



^EE / /A,,„G(Zi(a),Z2(a))i^y.,F,. (.,^:i^,^^,Z,(. 

Viip) 



(5.41) 



(5.42) 



-^liApi) = IKpiT' 

api 



-j—kApi) - (^i(pi)i ■j—h{pi))R<ikApi) 

api dpi 



(5.43) 



(2) (2) (2) (2) (2) 

We deal with J-j; 2 first. It may be split as J{ 2 2 1 + "^1 22 + "^i 2 3' where 



7(2) ._ 
'^1,2,1 •" 



X Cl 



^EE / /A,,.G(Z,(a),Z2(.))i?,.,,F,„ U^^,^^M-)) 
p) 



(5.44) 
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u s 1 



(2) ._ 
1,2,2 •- 



Li{a,s) L2{cr,s) 



,lj{a) 



and 



U S p 



J 



(2) ._ 
1,2,3 •" 



X Ci 



Li{a,s)\ d 



Li{a,s) L2{a,s) , 



-h,/^{pi) dsdpdpi. 



6 J dpi 

By virtue of Lemma |5.4| , ( |5.32| ) and the definition ( |5.15| ) , we obtain easily 

\nJ?l2C]\ = O(5^/'"'^0l|VG||^..((iR.)2)(n - t)EC, as 5 ^ 0. 



(5.45) 



The same argument also shows that |E[J-[2 3C]| is of the order of magnitude of the right hand side 
of (|5l5|) . 

Using Lemma 5.1 and the definition ( ^.151) we conclude that 



Li(a,s) L2{a,s) 



6 



,h{a)] =0(6). 



Therefore, |IE[J^^2 is equal, up to a term of order 0{5^ '^^){'^ ~ 0I|VG'IIl°°{(R'*)2)EC, to 



Li{(y, pY 



(5.46) 



t a 



X Z)y,,^Dy„„ci (^M^j ci ^^inrill^J |/.(a)|Z,^(a)C 
Let 6 < (2pi)i/(i-'>'i) and fix i. We may apply Lemma |5.3| , with 



ds dp. 



Z = D,^^^Gih{a),l2{cTm (^C^^^^^,=^^,Z.(^)J Ua)\k,p{a)C, 

Xi := Dy^^i3Dy.^aCi{x), X2 := ci(x), 

Li{a,s) Li{a,p) 
91 ■■= ^ , 92-= ^ , r = Ci3{p-a),n=Ci3{s-p), 

where C13 > depends only on N. We conclude that 



E 



1 



i,Q,,(3 



E 



e. ( a,^^^^,^^^^,h{a) ] A„aG(Zi(a),Z2(^)) (5.47) 



6 



t cr 



xUa)\k,f,{a)dlpR 



Li{a,s) - Li{a, p) 



5 

u s 



ds dp 



< C7Mri||VG||,.((^.p)E[C] 1 1 ^i/2^ gi3(^P^ ^) ^ ^1/2 

i cr 



Cl3{s-p) 
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ds dp. 
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Here we used the fact that 



Li{a,s) L2{a,s) 
1 ^' ^ ' ^ , t-t(^) 



/ (pi) Li{a,s) L2{a,s) \ 



The right hand side of ( ^.471) is of the form Ci5{5){u — t)\\V G\\ ^oo (^(^d-^2-^MC, where Ci5{5) vanishes, 
as 5 — > 0. The second term appearing on the left hand side equals to 



u 

i,a,f3i ^ ^ 



Li{a,s) L2i(7,s) 



li{a) A„aG(il(cT),Z2(cT)) 



(5.48) 



dp 



s-p 



li{a) dp 



\li{a)\C } ds 



i,a,l3 i 



Lijcr, s) L2{a,s) 
6 ' 6 ' 



li{a) ] Di^,^G{hia),l2{cT))daR ( ) Ua)\C 



s — a 



ds 



thanks to the fact that \7yR{0) = 0. The term appearing on the right hand side of ( |5.48| ) vanishes as 

(2) 

5 — > and, in consequence we have shown that |E[Jf2'iC]| = Ci6{6){u - t)\\VG\\ ^.oci^^jjdpjEC, where 
Ciq{6) vanishes, as (5 — > 0. 

We now estimate given by ( |5.41| ). Note that according to ( |5.43| ) and ( ^.16| ) we have 



t(2) _ 7(2) , j(2) 
•^1,1 — "^1,1,1 "^1,1,2' 



where 



with 



while 



u s 

X (s - pi)Ti^p ( pi, — - — , — - — ,Zi(cr) I dsdpi, 
Ti 2/1,^2,1) := |ir^ Fi{p,yi,y2,l) - (l, Fi {p,yi,y2,l)^ ^^l 

u s pi 

•^^2 ^=^EE/ / /A,,aG(Zi(a),Z2(a))Dy,,^F,,« (s. 



Li{a,s) L2{a,s) 



5 



viipi) y2{pi) 



li{p2) dsdpi dp2. 



(5.49) 



5 ' 5 

A straightforward computation, using Lemma ^ (note that -^^i,f3 ~ in ( 5.49D ), shows that 



the terms E[J{ 2 2C] ^-^icl E[J-[;^2 3C]) yields that E[J-[;^/ is equal, up to a term of the order Gn{6)(^ 



r(2) 



An application of Lemma pA, in the same fashion as it was done in the calculations concerning 



(2) 
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*)l|VG||ioo((Kd)2)E[C], where f^^Cn{6) = 0, to 



i,a,l3 



A,,aCr(ti(o-),t2(o-))-C'y,,/3-f^j,Q s, 7 , 7 



(5.50) 



t cr 



X J- i,/3 Pi, 7 , 7 , ) C 



ds dpi 



We denote 



\ 7:9 7 / 



Applying Lemmas |0| and 5^, as in (|5.46|) ( 5.47 ), we conclude that (|5.50| ) is equal, up to a term of 



r,cx)rriijd^2)E[<^], where ViinCis{S) = 0, to 



order Ci8(<5)(n-t)||VG||io 

u s 

^Yl j /(s-pi)IE{A„aG(Zi(a),Z2(a))A,(a,P,;Z(-))yi,«(P0C} dsfipl, (5.51) 

' t cr 

with Aj defined by ( |5.13| ), and Pi = (Zi(o", s), ^2(0', s), Zi((T, pi), -L2((T, pi), Zj((T)) . Note, however, 

(s, i:i(cr, s), Z2('7, s)) = E (^S^^'\ ii(cr, s), i:2(cr, sfj 



that for. G[47),4^i 



and 



; (a, Li{a, s),L2{a, s)) - E {a, (<t), y2(^)) I <CJ2 «) " ?/p(^) | < C7(s - a) < CJ^" 

p=i 



■71 



A similar estimate holds also for the terms containing Li{a, pi) and we conclude that the expression 



in ( 5.51 ) is equal, up to a term of order Ciq{5){u — t)||VG||^tx;((iRd)2)E[C], where lim Ci9((5) = 0, to 



U I s 

/ ^lD,^,aG{h{a),l2{a))e-ia) J {s - pi)V^,^{P,) dpi 



C > ds, (5.52) 



with 



G,(ct) ■.= e,ia,yi{a),y2ia),liay,l{-)). 
Note that, for s > t + 6^""^^ we have 

s s 



(5.53) 



E^',7,.^ (^T^'^^^O ^;.('^)^;7(^) -E^l7,7« ( 



S - Pi 



liicr) 



7.9 

Ei'^(^)i / PI 



dpi 



E ^",7,9-^ ^»(^)) 49(f^)47('^) - E ^",7,7^ (Pl '"»(^) 



7.9 



dpi 
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Using the fact that 



we obtain, upon the integration by parts performed in the first term on the utmost right hand side 
of ( 5.53 ), that this term equals to 



(r^i k,^{a) - J d^^R {pi k{a)) h,^ 



[a] dpi 



(5.54) 



pidl,y,yR(pili{a)) dpi 



5-71 



= J2 [^"^'^',7^ {^~^" U^)) - daR (S""'' Ua 

We have used here the fact that V-R(O) = and 



pidl .^ .^R{pi li{a) ) dpi 



7 

Summarizing the work done in this section, we have shown that 

|IE[/('^C]| < C2o(n-t)||VG|L=.((K.p)EC, 
where the constant C20 does not depend on 5 and G. 



(5.55) 



5.3.2 The terms E[/(2)C] and 

The calculations concerning these terms essentially follow the respective steps performed in the 
previous section so we only highlight their main points. First, we note that because li{p) — ~ 
,51/2-71 we have that E[/(2)(^] ig, up to a term C2i(5)('u-i)l|VG||ioo((Kd)2)E[C], where lirn C2i(5) = 0, 

equal to 



<5— 



Dy^^^G{h{a)M{<y))D,^,pFj,a ( 5,^,^,i,(a) 



E 

viip) y2{p) 



(5.56) 



ds dp. 
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Replacing p by a as the argument of hi') in (|5.56| ) needs a correction that is of order of 

magnitude 0{5^/^~^'^^){u-t)\\VG\\Loo((^^d)2)E[C], since 71 E (0, 1/8]. Next we note that (|]5|) equals 



E 



t a 



-y 

/ Li{a, p) L2{a, p) , . . , ^ 
X ( ' J ''i(^) I C 

u s 1 

E 



(5.57) 



ds dp 



A,,aG(Zi(p),Z2(p))i:'y,,T,A,,/3i^j,Q S, 7 , 7 ,Zj(cT, 



i,7 i,a,/3 t ^ 

Li(o-,p) L2{a,p) 



<52 

X ( P 



' 6 

u s 1 

E 



ds dp dv 



i,7 j,a,l3 J ^ 

i?i(u,cj,p) R2{v,a,p) 



Di^^aG{li{p),l2{p))Dy,^^Di^,l3Fj,a ( S, — ^ , , Zj (cr) 



ds dpdv. 



A simple argument using Lemma |5.4| , ( ^.3lD and (|5.32D shows that the second and third terms of 
(|53^ ) are both of order of magnitude 0{6^/'^~^'^^){u - i)||VG||ioo((Kd)2)E[C]. 



The first term, on the other hand, can be handled with the help of Lemma 5.3 in the same 
fashion as we have dealt with the term jj^2 i' given by ( 5.44| ) of Section 5.3.1 , and we obtain that 



|E[/(2)C]| < C22(<5)(n - t)||VG||^^((K.).)E[C], 



(5.58) 



where limC22('5) = 0. 

5^0 



Finally, concerning the limit of E[/'^'^^C] we note that by Lemma |5.4| we have 
ni^'k] < C23i6){u - t)||VG||^^((K.)2)E[C] + 5] J,,,- 



(5.59) 



where lim C23(5) = and 
5^0 



u s 

^^'1 '■=\T.I I E|A„;3A,,aG(Zi(a),Z2(a)) 



X Fj- Q s 



Li{a,s) L2{a,s) 



Li{a,p) L2{a,p) 



,li{cr) I C \ dsdp, 



First, let i / j and 25^-'^i Af < (2g)-^ Suppose also that s > t^f\ We have then 



\Li{a,s)-Lj{a,p)\ >--2M{s-a) > — 

q 2q 



on the event (with fixed a, (3) 

Li((j, s) L2{a,s) 



,lj{a) ] &i(s, , ,li{a) 
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xDu,(3Di^,aG{liia),l2ia))\ljia)\Ua)\ + 0. 
When, on the other hand, s,p G [0,t^^^], then we conclude from ( |5.3| ) that 

\Li{a,s) - L.j(a,p)\ > 7^5 > 7^(5 - a). 
Therefore llij] can be then estimated via Lemma ^.31 and Lemma by 



E[C]. 



(5.60) 



It obviously vanishes, as (5 — > 0. The second term in ( ^.601 ) arises from the contribution of s < 

When i = j we can use Lemma 5.3 in order to obtain 

< C2a{u - t)||v2G||ioo((Kd)2)E[C]. 

Summarizing, we conclude that 

|E[l(3)C]| < C25{n - t)||v2G||ioo((K.)2)E[C], 

where C25 can be chosen independently of 6 and G. Hence we conclude ( ^.29 ) and tightness follows. 

Suppose now that Q is any limiting measure of Q^^"-^ for a certain sequence (5^ — > 0, as n ^ +00. 
Coming back to ( 5.52 ) we conclude, using calculation ( ^.53D -( ^.54D , that the limit, as (5 ^ 0, of the 
expression on the left hand side of ( 5.52| ) equals to 



a^\s)Di^,^GiK^is),K2{s))\K,{s)\A,{s)C{K)ds 



Q{dK), 



(5.61) 



where 



Ai{s) ■.= Ai l^s,Yi{s),Y2is),Ki{s),K2{s)-Ki-),K{-) 

s 

Yi{s) := Xi + j ki{p) dp, i=l,2, 


+00 



Similarly, we calculate the limit, as (5 ^ 0, of E[/(^)C]. We know that only the limits of the terms 
Tj^j contribute. A straightforward computation shows that 



u 

lim^J,, = I I c%{s)Dy^,^D,^^pG{Kr{s),K2{s))H^'^{s)aK)ds 



Q{dK), 



where 



cS/3(s) ■■= -\Ki{s)\' J dlf,R{pKi{s))dp, 


H^'^is) ■.= eUs,Yi{s),Y2{s),K,{s)). 
Summarizing, we have shown that any limiting measure Q satisfies (|5.30| ) 
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5.4 The removal of cut-offs and the proof of weak convergence of {kf\-),k2\-)) 
Let 0ki,k2 •= ^ki ® 0k2 be the law of two independent copies of the diffusion given by (|4.1|) over 



C2{ki,k2) starting respectively at ki and k2. For a fixed M let 0^^^ be the law over Ci of any 
diffusion starting at a given ki G with the generator given by 

£(^^)F(k) = 5]a(^^)(k)a2^,, F(k) + j;6f)(k)afe,F(k), F G C^{R'^). 
p,q p 

Here ap^P{-), bp^^\-) are bounded and twice continuously differentiable, ap*g'*(k) = |kpZ)p^g(k), 
6j,*'^)(k) = |k|£;p(k) for M-^ < |k| < M. By virtue of Theorems 5.2.3 and 5.3.2 of HI we conclude 
that nj^f ^ is the unique probability measure such that 



z 

F{K{t)) - F(ki) - j 6^^^F{K{s))ds, t > 



is an (^M5L'*)^^^-martingale for any F G Cl{W^). We define Q^^^]^^ := o!^^ ^ ^ 

Let us briefly describe the strategy of the proof of weak convergence of (•))• First, 

for any ii' G C2 we define a certain stopping time W{K; M, N,p, q), see ( ^.63|) . The crucial property 
of that time is that the dynamics given by ([4.8D agrees with the dynamics of the truncated system 



(|5.16| ) up to W{-;M,N,p,q). We also show that any limiting measure Q{-;M,N,p,q) satisfies, up 



to the stopping time, the martingale problem associated with the diffusion given by l}ki,k2- This 
property allows to identify Q{-; M, N,p, q) with i3ki,k2 on the cr-algebra A^g'^ corresponding to the 



stopping time. The final step is to show that for sufficiently large A^, so that (5^) is satisfied, and 
sufficiently large M, as in ( |5.lD , the stopping time W{-; M, N,p, q) converges to infinity in nki,k2 as 
q — > +00 and p — > +00 (in that order), see (5.64). The weak convergence statement is a consequence 



of this property of the stopping time and it is shown in the calculation following ( ^.80]) . 



We introduce the following (A^2'*)t>o^stopping times. As before, for any K = {Ki,K2) such 
that K{t)j^O for alH > we define 

t 

Yj{t) := [ kj{s)ds. (5.62) 







Ap) Ap) 

''k ' ''k+1 



For such a K we let S{N,p) := lim Sn{N,p), where 

nf +00 

SniN,p) := inf |t > : for some A; > we have t G 

and KAtf) ■ KAt) < 1 h -, for some i G {1, 2} or j G {fc - 1, k] 

■' A* + 1 n 

If K is such that it becomes for some t we adopt of the convention that S{N,p) = +00. We let 
further T{M) := lim Tn{M), where 

nf +00 

r„(M) := inf > : \Ki{t)\ < ^ + ^, for some i G {1, 2}, 

or \Ki(t)\ > M - -, for some i G {1,2} 
n 
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Finally, for any Ri,R2 > and K = {Ki,K2) G C2{Ri,R2) we let U{p,q;K) := lim Un{p,q; K), 

nf +00 

V{p,q;K) := lim K), where 

nt+00 



Unip, q; K) := inf |t > : for some fc > 1, i G {1, 2} we have u G [0, 4^i], * G 



such that |yj(f) - Yi{u)\ < - + - I , 

q n J 

Vn{p,q;K):=\ni\t>-: inf |yi(t) - ^2(^)1 < - + -, or inf \Y2{t) - Yi{u)\ < - + - 
IP 0<M<t q n Q<u<t q n 

We adopt the convention that any of the above defined stopping times is infinite if the respective 
set of times over which it is determined is empty. 

Suppose that Tq > is an arbitrary deterministic time. Let 

W{M, N, p, q) := S{N, p) A T{M) A q) A V{p, q) and (5.63) 

B{M,N,p,q) := {S{N,p) AU{p,q) AV{p,q) <T{M) ATo}. 

We have B G A^g '^. According to Theorem 6.1.2 of ||2^ the measures i3ki,k2) ^k^/kj' Q(''^ -^'P' l) 
agree, when restricted to A^2'^. In what follows we show that 



lim lim Ok, k, [W{M, N, p, q) < +00] = 0. 

p— >+oo 5— >+oo ' 



The condition 



To < W{kf\.),k',">{.);M,N,p,q) = W{ir{-)A"'i-y,M,N,p,q) 



(S), 



(5.64) 



(5.65) 



implies ikf\s),k^2\s)) = (zS^^(s),z5'^^(s)) for s G [0,To]. We will use both ( p4| ) and (|]6|) to 
establish weak convergence of the laws of {kf \-),k2 \-)) over C([0,To]; ''td"^''^^ 
We start with the following simple observation. 



Lemma 5.6 With the choice of M as in (|5.l| ) we have 

nki,k.[r(M) = +cx)] = i. 



Proof. A simple calculation using Ito formula and Remark 4.1 shows that d|kj(t)| = 0, j = 1, 2 
which proves the lemma. □ 



Lemma 5.7 Under the assumptions of Theorem |4.4| we have 

lim U{p,q) = +oo,\/p, i3ki,k2 -a.s. 

(J— > + co 



(5.66) 



Proof. The proof is essentially the repetition of the argument from p^ ] pp. 60-61 so we only 
highlight its main points. It suffices only to show that 



lim U^'Hp,q) = +oo,\/p, Ok, 

(J— > + co 



a.s., 



(5.67) 



where C/(*)(p, g) := lim Un\p,q), 

nf +00 



Uli^p, q) := inf |t > : for some k > 1, we have u G [0, 4^2^], t G 
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such that 



Ki{s)ds 



1 1 
< - + - 
q n 



However, ( p.67| ) can be proved with the help of the argument contained in pp. 60-61 |15] so we omit 
the details here. We obtain from ( |5.67 ) 

lim U{p,q) = +oo,Vp, Hki.ka - a.s. 

However U{p,q) = U^^\p,q) A U^'^\p,q) and (5.66) follows. □ 
Let us denote by 



2) 



U) ._ 



U ^^(^ 



(5.f 



0<s<t 



and by 5^(2)^^) := [x : dist(x, 2)1^0 < r] the sausage, up to time t, of diameter r > around 
trajectory Yj{-). 

The next lemma shows that S{N,p) becomes infinite as p — > oo for each A^. 



Lemma 5.8 We have 



lim S(N,p) = -|-C!0,V A'', Qk ko ~ Q^-*- 

p— >+oo ' 



(5.69) 



Proof. The conclusion of the lemma is a consequence of the uniform continuity of paths of the 
diffusion on any finite time interval [0,T], which implies that 

lim min min kJt) • KJ&^) = 1, j = 1,2. 
P-+-4'')e[o,T]te[4-)4!)j 

□ 

Our next lemma shows that V{p, q) becomes infinite with p,q ^ +oo. 

Lemma 5.9 Suppose that N is as in (^^) and Ti,ri > are arbitrary. Then, one can find po,qo 
such that 

Q^,,^,,[S{N,p)AV{p,q)<Ti]<r], Vp > po, <? > go- (5.70) 

In order to prove this lemma we will need an auxiliary property of {Ki{-),Yi{-)). Let ki = |ki|. 
Note that the process {Ki{-),Yi{-)) is a diffusion on R'' x M"^, actually supported on ^^"^(O) x M'^, 
over (Ti , £2ki ) • Its generator is given by 

AA(k,x) := £F(k,x) + k- VxF(k,x). 

We denote by P(t,k, x; •) its transition probability. It satisfies the Fokker-Planck equation 

+ 00 

/ / /(^* = °' G C^((0,+cx3) X R"^ X R'^). (5.71) 

Rd ]Rd 

Lemma 5.10 Let t > 0, (k, x) G 5'^""'^ x R'^ (k = |k|^. Then, P(t, k, x;-) is absolutely continuous 
with respect to the Lebesgue measure on 5"^ ""^ x R'^, with the transition probability density p{t, k, x, •, •) 
that is a C°° -function. In particular, for any T, K,ri > there exists a constant C > such that 

max max_ P(t, k, x; 5"^"^ x yl) < C| A| (5.72) 



for any A C Bk{0) and A G B{R'^). 
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Proof. Let k := |k| and S^^^ : -Bf"^(0) ^ S"^^^ be given by 



Srd) := (/i,--- ,±\/fc^^,--- 1 = (^1, • • • , Id^i) G Sf-'(O), / := |1|. 

j-th component 

Define the measure P,^^\t,B x A) := k, x; x A), ^ G B{R'^), B G ^(5f^^(0)). The 

conclusion of the lemma holds if we can show that each measure Pi possesses a C°° smooth density. 
In what follows we consider only the case i = d and denote 5 := S^^\ Ps := Pj"^''. Note that the 
respective measure satisfies the equation 

{dt-N*)Ps{t,-)=Q 

in the distribution sense, with 

d-l d-l 

AAF(l,x) :=/c2 D,,,{\)dliF{\,^) + kY,E,m,n^,^) 

p,q=l p=l 



p=i ' ^ 



where Dp^q{l) = Dp^q{k^H, k'^Vk^ - P), Ep{l) = Ep{k~H, k^^Vk^ - P). It suffices therefore to 
prove that dt — M* is hypoelliptic in order to prove the Lemma. Note that 

d 

{dt -Ar*)F = Y, + XoF + a(l)F, VF G C^iB'jr^) 
p=i 



where 



d-l 

Xp{l) :=kY,Dl/^m„ P=lr-- ,d-l, 

q=l 



q=l ' ^ ^ q=l 

and a(-), ap(-), p = 1, • • • , d — 1 are C°°-smooth functions. It suffices therefore to prove that for any 
(i, I, x) G M X S'^~^ X M*^ the linear subspace £i,i,x of the tangent space to M X X M'^, spanned 
by the vector fields belonging to the Lie algebra £ generated by Xq, Xi, - ■ ■ , X^, is of dimension 2d. 
The (d - 1) X (d - 1) matrix D(-) := [I)p,q(l)], as weh as D^/2(-), is non-degenerate in B'^, ^(0) due 
to Proposition (actually it degenerates in the limit when I approaches 5-B^~^(0)). Hence the 
vectors di^ G i2t,i,x, p = 1, • • • ,d-l. 
Note also that 
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where bp{-), p = 1, ■ ■ ■ ,d — 1 are C°^-smooth functions. Hence, d^^ + lq{k'^ — ^^'^dx^ G '2t,i,x, 



g = 1, • • • , d — 1. Furthermore, 



p=l 



trD(l) + (D(l)l,l)i,.(A;^-/^) 



2\-l 



9=1 



s d-1 



where Cp(-),dp(-), p = I,-- - ,d — 1 are C°°-smooth functions. We can conclude therefore that 
G •Ctj^xi hence also dx^ € -C^^i^xi p = 1, • • • , d — 1 and finally we also get dt G £f,i,x5 so that the 

proof of Lemma 5.10 is complete. □ 
Proof of Lemma |5.9| . Let 

A{N,v) := [S{N,p)>Ti + l]. 



Choose p sufficiently large so that 



£1w,m[MN,p)] > l-r?/2. 



This can be done thanks to the continuity property of diffusion paths. For any (A'i(-), i^2(")) ^ 
A{N,p) we have 



Yi 



Y2{p) 



> Xn{p), and 



Yiip) 



> A7v(p) 



(5.73) 



for all p G [0, according to (|5.4|) (see (|5.62 ) for the definition of Yi{-), i = 1,2 ). Recall also 
that ^f\Ki), i = 1,2 are defined by (jS^e^ ). 

Let F(i)(p,g;K) := lim yi^^(p, g; i^), where 

y«(p,(?;if) := inf |t > - : dist(yi(t), 2)f ) < - + -] 
IP q n ) 

and likewise we introduce V^'^\p,q;K) := lim Vn'^\q;K), with 

n— ^+00 

K(')(p,(Z;i^) := inf > - : dist(y2(t), 2)1'^) < - + -| 
[ p q nj 

Note that V{p,q;K) := T/(i)(p, g; i^) A ^. 

The conclusion of Lemma |5.9| is then a consequence of the following. 



(5.74) 



Lemma 5.11 For any N sufficiently large so that ( |5.3D holds and p > 1 we have 

lim V^^\p,q;K) = +0O Oki ka - on ^(iV,p), Vki, k2 / 0, i = 1, 2. 

Proof. With no loss of generality we assume that i = 1. Note that obviously 
V'^^\p,q';K) > V^^\p,q;K) for q' > q and ah K G C2{ki,k2). 
For any K2 we denote by 

A{N,p-K2) ■■= [Ki : {Ki,K2) G A{N,p)]. 
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It suffices to show that for £)k„-a.s. K2 we have 



hm V'^^\p,q;Ki,K2) = +00, 0ki - a.s. on A{N,p-K2). 

9— >+CXD 



Let us denote 



B(t,x;K2) := [Ki : ^(t; x) - YsCp)! > Ajv(p), P G [0, 1/p] ] . 
Note that A{N,p; K2) Q B (1,0; K2) , according to (|5J^ ). 



0, ~ ill -^2 



Let Ti > be arbitrary. We show that 



hm Oki 

g— »+oo 



(5.75) 



The expression under the hmit in (5.75) can be estimated by 



inf distiYi{u),^^^^)<-,B(-,0;K2) 
ne[o,Ti] '1 q \P J 



St;^>cll/P>\^\>^r^{p)] 



P [ -,ki,0,dk,dx ) Hk 



inf 



if dist (yi(n;x),2)5?)) < B{0,^;K-. 



ne[0,Ti-l/p] 



Here we used the Markov property of the process (i^i, Yi). ( |5.75 ) fohows if we can show that 



hm £2k 

(J— > + oo 



inf dist (yi(n:x),2}^^h < - 



(5.76) 



for every k S ^ and x satisfying 1/p > |x| > Xn{p), £2k2 ~ ^'•S- ^2- 
Suppose first that r/i := ^dist(x, S)^^"*) > 0. Then, 



£2k 



0<n<7?i 



inf dist (yi(n;x),2)^^^ ) > - 



1, Vg>4r/-^ 



(5.77) 



Note that the expression under the hmit on the left hand side of (^.76|) can be estimated by 



Hk 



iiif dist (Fi(jyg;x), 2)52^) 



2 

< - 



(5.78) 



<{Ti + l)q max nk[>l(j7g;x) Gi?2/g(2)Ti^) 



The right hand side of ( |5.78D can be estimated, with the help of Lemma [5.10| , by 

(recall that l2(-) is of C^-class, with |y2(")l ^ 1) ^^^^ ( 5.75 ) follows, provided we can prove that 



(2) 



dist (x,2)fM = 







(5.79) 



for 1/p > |x| > XnIp) Recall that |x — l2(p)| > '^Nip), P G [0, For any p > we can estimate 
therefore the left hand side of ( 5.79| ) by 

Ti + 1 



4-1 



max £lk2 [ \Y2{jp) - x| < 2p] < C{p){Ti + l)p' 

P l/p<jp<Ti 

for some constant C(p) > depending only on p. Since the last inequality holds for all p > we 
conclude (|]7|). □ 

An immediate consequence of Lemmas 5^, 5^ and 5^ is the following. 
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Corollary 5.12 For any M,e > there exist find sufficiently large p,q and N so that 

Qki,kj5(M,A^,p,g)] < e. 

T 

Choose any Tq > and F a bounded and continuous functional over C2 that is M2 "-measurable. 
We show that 



limsupE F{k[^\-),k'°'i-)) < / F(K(-))Qki,k2(t^i^)- 



(5.80) 



This, in fact, implies weak convergence of the laws of {k[^\-) , k^2^\-)) over C2, as (5 — > 0. 

Fix 77 > and choose M > such that M — 1 satisfies (^.11). Then, by virtue of Lemma |5.6| 



0ki,k2[T(M-l) <ro] =0. 



(5.81) 
(5.82) 



Let p, q be such that 

£lkiM2lB{M,N,p,q)] < 77. 

Note that B{M - 1,N - l,p,q - 1) O B{M, N,p, q). 

Let 6n — > 0, then we can choose a subsequence, that we still denote as {6n), such that the laws 
of converge over C2, as n ^ +00, to a certain Q{-; M, N,p, q). We have 



< limsupE 

n— >+oo 



limsupE F{kf"\-),k^^"\-)) 

n— >+oo L 

j{Sn)( \ j{5n). 



(5.83) 



F{l\°-\-)a'2'''{-)y,W{lf-\-),l^^-\-y,M - 1,N - l,p,q - 1) > To 
Fik^^-\.),k^^"\-))-Wilf-\-)/^-\-y,M-l,N-l,p,q-l)<T, 



+ lim sup 



E 



The second term on the right hand side of ( p. 83 ) can be estimated by 

flimsupF[r(zP(-),Z^^(-);M-l) <roj +nki,k2 [B{M-l,N-l,p,q-l)]) (5.84) 



Note also that 



lim sup 

n— >+oo 



T{l[^"\-),l[^"\-y,M -1) >To -P VI/(Zp^(-),4''"^(-);Af - l,Af- 1) >To 



{5n) / \ l(5„) , 



< lim sup P 

n— ++00 



(5.85) 



(5.86) 



e BiM -1,N -l,p,q-l) 
<£}ki,k2 [:B(A/-l,iV-l,p,g-l)] 
and hence 

limsupP \T{lf''\-),l^^"\-);M - 1) < Tq 

n— >+oo L 

= 1 - hminf P \T{lf-\-),l^^-\-yM - 1) > To] < Ok^.k^ [W{K;M -1,N- l,p,q- 1) < Tq] + r/. 

n— >+oo L J 

The first expression on the utmost right hand side of ( ^.861 ) is less than or equal to 

Oki,k. [T(i^;M-l) <ro]+£3ki,k2 [B{M -1,N -l,p,q-l)] < r/ (5.87) 
according to (5.81) and ( |5.82| ). Summarizing, the expression in ( |5.84| ) can be estimated by 27/11^112,00 
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The first term on the right hand side of ( |5.83| ) can be estimated by 



+ \\F\\L^£l^,,^,[W{K;M,N,p,q)<To] < J F{K{-))Q^,^^,{dK) + 2t]\\F\\l^ . 

The last estimate follows from an analogous estimate to ( ^.871) . Summarizing, since ry > is arbitrary 
we conclude (|5.8C1|). 



A Proof of Lemma |273 

We define 



4(^,xo)=vf'^-vP = |, 



,ik-(4-y) 



and split dg = Yl'j=i '^e according to the decomposition 



+ ([/^(xo + ^^^^,k) - C/'^(xo,k)) + ([/^(xo,k) - H^,^(xo,k)). 



roSo(y)H (A.1) 



,k)) 



Here = Y. UgUg is the semi-classical approximation of , with the solutions of the Liouville 
equations. The last term may be estimated as 



df(^,xo)fdxo < 



e^'^-^[[/'^(xo,k) - VF/(xo,k)]roSo(k) 



dk 



dxo 



< C\\So\\l2 j ||C/^(xo,k)-T^,^(xo,k)||2dxodk^O 
as fC^ 3 (e, 5) — > with C independent of ^. The Fourier transform of the first term in xq is 



-,k)roSo(y)- 



2 

I e^i^-(^-^)+-P-(«+^)/2/,^(p,k)roSo(y)|;^ = I e'('^+^P/2)-V.^(p,k)roSo(k 



(27r)' 



ep dk 



where /e(x,k) = Wj^(x,k) - C/'^(x,k). Therefore we have using the Cauchy-Schwartz inequality 

I \\4'\tp)fdp < C||Sok2||/,^||^2 ^ 

5 2 

as /C^j 3 (e, 6) ^ with C independent of ^. Finally, the Fourier transform of de is 

d'Atp) = /e-P--+^''-(«--)[C/^(xo + ^^,k)-[/^(xo,k)]roSo(y)^^^ 



e*ep-{4+y)/2 _ I 



C/^(p,k)roSo(y) 



e^''<U\p,k)To [e-P-«/2So(k - ^) - So(k) 



dydk 
dk 



(27r)^' 
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We write e*^P'^/^ = (e*^P'^/^ - 1) + 1 and decompose (iP(^;p) as /i(^;p) +-^2(^;p) accordingly. We 
have for the second term 

/ \\l2{^;p)fdp<C I (I ||;7^(p,k)||||So(k-^)-So(k)||^ydp 
<c|(| \\u\p,k)fdk^ (/ ||So(l-^)-So(l)fdl) dp. 

Note that 

I ||So(l - ^) - So(l)|pdl < e^llplp I (^1^' ||vSo(l - eps)\\^ds^ dl < e^HpfUvSoHi^ 
and hence 

/ ||/2(^;p)fdp<e'||VSo||i.r^||^.-0 

as ]Cj_i B (e, (5) — > according to Lemma 3.6. It remains to bound the norm of /i(p; We derive 
two estimates according as ^ is small or large. The first estimate is 

|/i(^;p)l < Cy"e^i^-^|[/^(p,k)|e|p.^||So(k-^)|^ 
< Cemu' {pMlLlipMllSoh, 



so that 



I |/i(^;p)l'^ip < e2|^|2||c/^(p,k)||^,(^,^^,(^,))||So||i. 

At the same time using integrations by parts we get 

/i(^;p)= |^e^'^-^^.Vk(f7^(p,k)roSo(k-^))(e-P-^/2-l) 



This shows that 



J |ii(^;p)l><|^l|f/|li.(K.^HMM^))ll(i + l-l')'^'So||i. 



With these estimates, we obtain that 



J |/i(^;p)P^ip<Cmin(/if|^|M^| 



with /if ^ as /C^ 9 (e, 6) 0. This implies that / p)pdp 0, hence / p)\'^dp 

as /C^ 3 {e, 6) ^ uniformly with respect to ^ G M'^ and concludes the proof of Lemma I 



B Proof of Lemma 13.3 



We may recast (1.5) as 



Thanks to calculations of the form ( |3.10| ), this is equivalent to the equation 



£^ + P,^'^(x,eDjvf = 0, (B.2) 
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where the symbol is given by (|3.8D . We recall that the pseudo-differential Weyl operator 
P^(x, eD) associated to a symbol P(x, k) is defined by Weyl's quantization rule 



(B.3) 



The fact that ( [B.2D is equivalent to ( |B.1| ) is verified by a straightforward calculation: 

d^dy 



^e / e 



9 



X + y 



e / C5 



x + y 



y=x 



9vf e dcs (x) 



(27r)^ 

ec,(x)«..^ + |%Wo^vf(x) 



DWt{y){-^^5{^-y)]dy 



and now (|B.2| ) follows because P^'^(x)v^(x) = Pj^(x)v^(x) since i^i (x) is independent of k. 
The associated Cauchy problem for the Wigner transform with a fixed Q is given by 



+ W[Pt'^i^,eD^)vi,vi] + Ty[vf,P,''^(x,eDjvf] = 



W^/(0,x,k)) = W^,^(0,x,k;C), 
where the Wigner transform of two different fields is defined by 



(B.4) 



(B.5) 



We deduce from the definitions of We and P^ that 



W[P,^'^(x,eDJvf,vf](x,k) 
I e^'^-y(Pf (x,el).)vf)(x 



g*k-ygi(x-^-z)-qp5|^ 



X-^+Z 



2 



.q)vf(z)vf(x + 



ey dydzdq 



2 ' (27r)2'>! 



(B.6) 



(vrV2) 



n^(y,q)T^,^(z,p)e^^-^^^P^^^^ 



2d 



Moreover, the matrix 14^^^ is self-adjoint, while W^[/e, ge] = W*[g, f] for any pair of functions / and g, 
and the symbol Pg is skew-symmetric. Thus ( p.6D and ( p.4| ) imply that the pure Wigner transform 
We satisfies ( 3.11| ), and hence so does We. Moreover, the function (p satisfies an anti-symmetry 
relation 

(/)(x, z, k, p; y, q) = -0(z, x, p, k; y, q). 
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Then, using the fact that is self-adjoint we obtain 
Tv{{4w^)Wt)d^dk 

= I TV(P,^(y, C,)W!{., p)W^{^, k)e^f - W,{z, p)P^{y, ci)W^{., k)e-'^) '^^'^^dy f ^^Pd^ 

= / Tr{P^{y,ci)W^{z,p)W^{^,k)e'^ - P/(y, q)W^,^(z, p)I4f (x, k)e^^) = 0, 

where we interchanged x •s-^ z and k •f-^ p in the second term on the last line, and used the anti- 
symmetry of (f>. This implies conservation of the L -norm Note that follows immedi- 
ately from ( p.ll| ) and the proof of Lemma |3.3| is complete. 

C Regularity of the Liouville equations 

We prove Lemma in this Appendix. We recall that the functions satisfy the evolution equations 

^ + {A^,n^} = (C.l) 

= Uq{t = 0) = Tv[UqWoUq]. 

These equations can be solved by following the Hamiltonian flow generated by A^. More precisely, 
let us define for T, x, k given, the trajectories 

dX{t) dXl 



dt 5k 



(X(t),K(t)), X(0)=x 



^ = ^(X(t),K(t)), K(0)=k. (C.2) 

Then solution of (|C.l| ) is given by 

ul{T, X, k) = n°(X(r, X, k), K(T, X, k)). (C.3) 

The flow ( [4.5| ) preserves the Hamiltonian Aq(x, k) and the initial data is supported on a compact 
set S. Therefore the set 

S= \^ supp Ug(t, X, k) 
f>o,5e(o,i] 

is compact because the speed C5(x) is uniformly bounded from above and below. Furthermore 

Vul = D'*Vul \\u,{t)\\m < HWmWD'm^ 

where D^{t,x,'k) is the Jacobian matrix, Z)j'* = dZf/dzj, with det D^{t) = 1, Z = (X, K), and 
z = (x, k) . To simplify notation, we do not write explicitly the dependence of and its derivatives 
with respect to the eigenvalue label q in the sequel. Here we define 

\\D^\\oo= sup fT¥[I)^(x,k)L»^*(x,k)])'^'. 
(x,k)e-s ^ ^ 

More generally, given a tensor Tj-^j^^^j^ we denote 

1/2 



|r||oo= sup I >^ iTj 
(x,k)e5 
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We will also use the matrix norm |^| that is dual to the Euclidean norm on M.'^ and is equal to the 
square root of the largest eigenvalue of the matrix AA* , and denote 



Furthermore, we have 



so that 



l^loo = sup |A(x,k)|. 

(x,k)e5 



dzjdzp dzjdzp dzm 



+ 



dZj dzp dzmdzr 



dzjdzp 



^ .. d^Zj d^Zf 9ng ^ ^ dZjdZ^. g^n" dZf dZ^^ Q^n" 



dzjdzp dzm dzjdzp dzi dzj dzp dzmdzr dzj dzp dzidzg 



< 



dzjdzp 



dzidzs 



5 |2 
mj\ 



and hence 



00 5 



with D 



2di 



d'^Z^/dzjdzi. We observe that 



Q3„.5 



U" 



d^Zt 92^^ Qzf a^nO d^zt dzi d-'ul 



+ 



+ 



dzjdzpdzs dzjdzpdzg dzm dzjdzp dzs dz^dzi dzjdzg dzp dzmdzr 



+ 



dzi d-'zt d\ 



+ 



dzi dzt dzf d'^u'l 



dzj dzpdzs dzmdzr dzj dzp dzg dzmdzrdzi 



Therefore we have 

E 



dzjdzpdzg 



< 



m,j,p,s 



dzjdzpdzs 



5\\m\L E 

m,r,l 



dZmdZrdzi 



SO that 



+ 3H\\H^\\DlmL^\\D' 



+ H\\m\\D' 



with -DgJ^p = d^Z^/dzjdzidzp. 



It thus remains to estimate the matrices D^, D2, and D^. The matrix satisfies the differential 



equation 



dkidxj dkidk. 



d\o) =i. 



\ dxidxj dxidkj I 
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Therefore we have 

di 



so that \\D^{t)\\oo < exp {\F^\oot) and hence 



< HWm exp IF-^ loci 



Differentiating (^]^) once again we obtain 
so that along each characteristic 



where i^2ijfc ~ dF^j/dzk- Furthermore, initially at t = we have D2{0) = 0. Therefore we obtain 

II r)§f+\\\..^ <r 



p2WllLoc.<^^exp(2|F^Ut). 



and thus 



Similarly, the tensor satisfies the ordinary differential equation 

3,jkm jpS T^S,p T^S,n j^S,l . riS,n j^5,l . rnS r>S,n j^5,l . rpS nS,n t^5,1 . rpS t^S,1 
- i'z,ilnp^mJ->k + ^2,iln^2,kra^j + ^2,ilnJ^k ^2,jm + ^2,iin^m J->2Jk + ^il^^jkm 

SO that along each characteristic 

■'- II n'5||2 _ 775 n<5iP n'''' n*^'* _l t?^ n^'^ n^,!- n^A 

2dt - ^3,ilnp^m ^k ^Z,jkm^^2,iln^2,km j ^3,jkm 

"r-^2,iZn-^A: ^2,jm-^3,jkm "i" ^2,iln^m ^2,jk-^3,jkm "i" ^ il-^3,jkm^3,jkm 

< ||F|||||Z)^f 111)3^11 +3||F2^||p^||p2'llP^II + li^'lll^^^ 
where F^-j^^ = dF2 ^ji^/dzn, and at t = we have 1)3(0) = 0. Therefore we obtain 

?<5|| Q||J?5||2 



\Dm\oo < + exp(3|F^Ut) 



1^ Mjk + ^JSI^^ ||nO|| + 3^^||n0|| + IkOi 



and thus 

^"^qKi^JWHi ^ [yiF^^ ^ |^5|2^ J ll"gll/ri |F^|oo ir"glli/3 

This completes the proof of Lemma ^!6| because 75 = |-F'^|oo- 



exp(3|F^|oot) 
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